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^ (54) Title: METHOD AND SYSTEM FOR PREDICTING NUCLEIC ACID HYBRIDIZATION THERMODYNAMICS AND 
^ COMPUTER-READABLE STORAGE MEDIUM FOR USE THEREIN 

^2 (57) Abstract: Method and system to predict and optimize probe-target hybridization are provided. The method may be imple- 
mented using six interactive, interrelated, software modules. Module 1 predicts the hybridization thermodynamics of a duplex given 

* — ; the two strands. Module 2 finds the best primer of a given length binding to a given target. Module 3 executes a primer walk to find 
alternative binding sites of a given primer on a given target. Module 5 is a combination of Modules 2 and 3. Module 6 finds the 
alternative binding sites of a given primer on a given target (Module 3) and calculates the concentration of target with primer bound 
at primary and alternative sites. Module 7 is a combination of Modules 2 and 5 and also calculates the various concentrations. The 

^ six modules can be operated either through an interactive user interface or using batch file submission as provided by Module 4. The 

^ program is suited to predict DNA/DNA, RNA/RNA, and RNA/DNA systems. 
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METHOD AND SYSTEM FOR PREDICTING NUCLEIC 
ACID HYBRIDIZATION THERMODYNAMICS AND 
COMPUTER-READABLE STORAGE MEDIUM FOR USE THEREIN 

BACKGROUND OF THE INVENTION 

5 1 . Field of the Invention 

This invention relates to methods and systems for predicting nucleic 
acid hybridization thermodynamics and computer-readable storage medium for use 
therein. 

2. Background Art 

10 Improvement of the efficiency of hybridization-based techniques 

requires the optimization of the binding between two sequences. Accurate 
prediction of the thermodynamics allows optimal choice of the sequences, 
temperature, and salt conditions. Hence, the prediction of nucleic acid 
thermodynamics is important to optimize techniques like PCR (Saiki et al., 1988), 

15 Southern and Northern blotting (Southern, 1975), antigene targeting (Freier, 1993), 
and Kunkel site-directed mutagenesis (Kunkel et al., 1987). 
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Hybridization prediction is also important for designing DNA 
microchips that have a wide field of application ranging from diagnostics (Hacia, 
1999; Yershov et al., 1996) to gene expression analysis (Ferea et aL, 1999) and 
drug discovery (Debouk and Goodfellow, 1999). Microchips contain a large 
5 number of DNA probe sequences that have to be designed to specifically hybridize 
target sequences in a pool of DNA fragments. First, a DNA probe should be 
designed to bind to only one site of only one DNA target. Second, the different 
DNA probe sequences need to hybridize to their targets under the same temperature 
and solution conditions. Moreover, in sequencing by hybridization (Fodor et al., 
10 1993; Mirzabekov, 1994) where microchips are used to determine the sequence of 
given DNA, one has to be able to know hybridization thermodynamics to 
discriminate signals resulting from perfectly matched and mismatched probe/target 
hybridizations. 

Another widely used technique that requires hybridization prediction 
15 is the fluorescence in situ hybridization (FISH) technique (Gall and Pardue, 1969). 
In this technique, a fluorescent tagged nucleic acid probe is designed to specifically 
hybridize cellular or tissue section nucleic acids. The target of these probes can 
either be endogenous DNA, messenger RNA or viral and bacterial sequences. 

Therefore, FISH is used to monitor gene expression (McNicol and 
20 Farquharson, 1997), detect infectious agents (Bashir et al., 1994; McNicol and 
Farquharson, 1997; Pollanen et al., 1993), study cell cycle (McNicol and 
Farquharson, 1997), map chromosomes and study nuclear architecture (Heng et al., 
1997). It was also determined that a set of probes can be used simultaneously 
(multiFISH) to detect different loci (Pagon, 1997). Once again, prediction of 
25 hybridization is essential to insure specificity. Nucleic-acid hybridization prediction 
is also important for the design of oligonucleotide aptamers or antisense 
oligonucleotides (Cohen, 1992) that can be used for various therapeutic applications. 
A new type of probes known as molecular beacons (Bonnet et al., 1999; Tyagi et 
al. 1998) that are very specific has been developed and shown to be efficient for 
30 mutation analysis (Giensendorf et al., 1998) and multiplex detection of single 
nucleotide variations (Marras et al., 1999). The design and prediction of the 
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thermodynamics of these beacons is helped by hybridization thermodynamics 
prediction (Bonnet et al., 1999). Accurate prediction of hybridization is also 
important for the practical realization of DNA-based or more generally nucleic acid- 
based computers. (Adleman, L.M., 1994). 

5 The development of molecular biology techniques based on 

hybridization (PCR, FISH, DNA microchips, etc.) has resulted in a need for 
efficient automated ways to design probes and primers. In the last decade, 
numerous algorithms have been developed to optimize the design of primers and 
probes for various applications (Rychlik and Rhoads, 1989)(Breslauer et al., 1986; 

10 Chen and Zhu, 1997; Dopazo et al., 1993; Haas et al., 1998; Hillier and Green, 
1991; Hyndman et al., 1996; Li et al., 1997; Link et al., 1997; Pesole et al., 1998; 
Proutski and Holmes, 1996). Numerous unpublished software to predict primers 
are also made available by research groups and biotech companies on the World 
Wide Web (Primer3 from the Whitehead Institute for Biomedical Research. Primer 

15 Express™ from PE Biosystems, DNAstar from IDT, etc.) 

There are currently many software packages on the market for DNA 
primer design including: OLIGO, PRIMER PREMIER, OSP, GCG, PrimerMaster, 
and Primo. None of the current programs, however, were written by experts in 
DNA thermodynamics; thus, there are many improvements that can be made. 

20 Nearly all of the current software packages contain mistakes that result from a lack 
of understanding of the underlying theory of DNA hybridization. PCR is a fairly 
robust process and thus even crude programs make predictions that work 90-95% 
of the time. Multiplex PCR primer design, however, is not at all trivial and 
detailed knowledge of the physical chemistry of DNA hybridization, and the 

25 availability of an accurate thermodynamic database are essential to reliable design 
of multiplex PCR primers. In multiplex PCR, several primers must be designed to 
specifically bind to different sites on target DNA at a given temperature with 
minimal background binding to mismatch sites and with minimal cross- 
hybridizations between pairs of primers. The design of molecular beacons for DNA 

30 oligonucleotide arrays is also very challenging because of the complex competing 
equilibria. 
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Most of the existing programs aiming at finding an optimum probe 
that binds a specific location on a target, however, do not include accurate stability 
rules for hybridization and neglect or poorly approximate competitive binding sites, 
strand folding and strand dimerization. 

5 U.S. Patent Nos. 5,593,834 and 6,027,884 to Lane et al. disclose 

methods to design and construct DNA sequences with selected reaction attributes. 

In summary, prediction of nucleic acids thermodynamics is important 
to optimize various molecular biology techniques including multiplex PCR, DNA 
microchips, molecular beacons, and fluorescence in situ hybridization. Most of the 
10 available programs for probe design do not include a complete parameterization and 
often do not account for mismatches. Moreover, single strand folding is not taken 
into account, which often leads to inaccurate predictions. 

SUMMARY OF THE INVENTION 

An object of the invention is to provide a method and system for 
15 predicting nucleic acid hybridization thermodynamics and computer-readable storage 
medium for use therein wherein the invention utilizes a thermodynamically rigorous 
approach to evaluate the quality of probes and simulate probe/target hybridization. 

Another object of the invention is to provide a method and system for 
predicting nucleic acid hybridization thermodynamics and computer-readable storage 
20 medium for use therein wherein the invention also takes into account single strand 
folding thermodynamics to calculate effective hybridization thermodynamics. 

In carrying out the above objects and other objects of the present 
invention, a method for predicting nucleic acid hybridization thermodynamics is 
provided. The method includes providing a database of thermodynamic parameters, 
25 receiving hybridization information which represents at least one sequence, 
receiving correction data, receiving a first set of data which represents hybridization 
conditions, and calculating hybridization thermodynamics including net 
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hybridization thermodynamics based on the 'hybridization " information, the 
thermodynamic parameters, the correction data and the first set of data. 

The hybridization thermodynamics of individual single stranded, 
bimolecular and higher order complexes may be statistically weighted in a numerical 
5 process and the equilibrium concentration of each species is output. 

The correction data may include folding correction data and/or linear 
correction data. 

The thermodynamic parameters may include DNA thermodynamic 

parameters. 

10 The DNA thermodynamic parameters may include dangling end 

parameters and/or coaxial stacking parameters. 



The DNA thermodynamic parameters may further include terminal 
mismatch parameters. 

The thermodynamic parameters may include RNA thermodynamic 
15 parameters and/or hybrid DNA/RNA thermodynamic parameters. 

The thermodynamic parameters may further include DNA loop 
thermodynamic parameters. 



The hybridization information may represent top and bottom strand 
sequences which form a duplex and wherein the hybridization thermodynamics are 
20 calculated for the duplex. 

The hybridization information may further represent at least a section 
of a target and a length of at least one primer or probe complimentary to the target. 
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The hybridization thermodynamics may be calculated for a plurality 
of primers or probes complimentary to the target. 

The hybridization information may represents at least a section of a 
target and a primer or probe. 

5 A length of the target may be longer than a length of the primer or 

probe and wherein the hybridization thermodynamics are calculated for a best 
target/primer or target/probe complex and for competitive mismatch complexes. 

Hybridization information may represent at least a section of a target 
and a primer or probe and wherein a length of a target is longer than the length of 
10 the primer or probe and wherein the hybridization thermodynamics are calculated 
for a best target/primer or target/probe complex and for competitive target/primer 
or target/probe complexes. 

The method may further include calculating concentration of each 
species in a solution at a plurality of temperatures. 

15 Hybridization information may also represent a primer or probe and 

wherein the length of the target is longer than a length of the primer or probe and 
wherein the hybridization thermodynamics are calculated for a best target/primer 
or target/probe complex and for competitive mismatch complexes and wherein the 
method may further comprise calculating concentration of every species in a 

20 solution at a plurality of temperatures. 

The hybridization thermodynamics may be calculated for at least two 
best target/primer or target/probe complexes and for their corresponding 
competitive mismatch complexes and wherein the method may further comprise 
correcting for any interactions between the at least two best target/primer or 
25 target/probe complexes and their components. 
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Further in carrying out the above objects and other objects of the 
present invention, a system for predicting nucleic acid hybridization 
thermodynamics is provided. The system includes a database of thermodynamics 
parameters, means for receiving hybridization information which represents at least 
5 one sequence, and means for receiving correction data. The system further includes 
receiving a first set of data which represents hybridization conditions, and means 
for calculating hybridization thermodynamics including net hybridization 
thermodynamics based on the hybridization information, the thermodynamic 
parameters, the correction data and the first set of data. 

10 The hybridization thermodynamics of individual single stranded, 

bimolecular and higher order complexes may be statistically weighted in a numerical 
process and the equilibrium concentration of each species is output. 

The correction data may include folding correction data and/or linear 
correction data. 

15 The thermodynamic parameters may include DNA thermodynamic 

parameters such as dangling end parameters. 

The DNA thermodynamic parameters may include coaxial stacking 
parameters and/or terminal mismatch parameters. 

The thermodynamic parameters may include RNA thermodynamic 
20 parameters and/or hybrid DNA/RNA thermodynamic parameters. 

The thermodynamic parameters may further include DNA loop 
thermodynamic parameters. 

The hybridization information may represent top and bottom strand 
sequences which form a duplex and wherein the hybridization thermodynamics are 
25 calculated for the duplex. 
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The hybridization information may also represent at least a section 
of a target and a length of at least one primer or probe complimentary to the target. 

The hybridization thermodynamics may be calculated for a plurality 
of primers or probes complimentary to the target. 

5 The hybridization information may represent at least a section of a 

target and a primer or probe. 

A length of the target may be longer than a length of the primer or 
probe and wherein the hybridization thermodynamics are calculated for a best 
target/primer or target/probe complex and for competitive mismatch complexes. 

10 Hybridization information may represent at least a section of a target 

and a primer or probe and wherein a length of a target is longer than the length of 
the primer or probe and wherein the hybridization thermodynamics are calculated 
for a best target/primer or target/probe complex and for competitive target/primer 
or target/probe complexes. 

15 The system may further include means for calculating concentration 

of each species in a solution at a plurality of temperatures. 

Hybridization information may also represent a primer or probe and 
wherein the length of the target is longer than a length of the primer or probe and 
wherein the hybridization thermodynamics are calculated for a best target/primer 
20 or target/probe complex and for competitive mismatch complexes and wherein the 
system may further comprise means for calculating concentration of every species 
in a solution at a plurality of temperatures. 

The hybridization thermodynamics may be calculated for at least two 
best target/primer or target/probe complexes and- for their corresponding 
25 competitive mismatch complexes and wherein the system may further comprise 
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means for correcting for any interactions between the at least two best target/primer 
or target/probe complexes and their components. 

Still further in carrying out the above objects and other objects of the 
present invention, a computer-readable storage medium having stored therein a 
5 database of thermodynamics parameters and a computer program are provided. The 
computer program executes the steps of; a) receiving hybridization information 
which represents at least one sequence; b) receiving correction data; c) receiving a 
first set of data which represents hybridization conditions; and d) calculating 
hybridization thermodynamics based including net hybridization thermodynamics 
10 based on the hybridization information, the thermodynamic parameters, the 
correction data and the first set of data. 

The hybridization thermodynamics of individual single stranded, 
bimolecular and higher order complexes may be statistically weighted in a numerical 
process and the equilibrium concentration of each species is output. 

15 , The correction data may include folding correction data and/or linear 

correction data. 

The thermodynamic parameters may include DNA thermodynamic 

parameters. 

The DNA thermodynamic parameters may include dangling end 
20 parameters and/or coaxial stacking parameters. 

The DNA thermodynamic parameters may further include terminal 
mismatch parameters. 

The thermodynamic parameters may include RNA thermodynamic 
parameters and/or hybrid DNA/RNA thermodynamic parameters. 
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The thermodynamic parameters may further include DNA loop 
thermodynamic parameters. 

The hybridization information may represent top and bottom strand 
sequences which form a duplex and wherein the hybridization thermodynamics are 
5 calculated for the duplex. 

The hybridization information may represent at least a section of a 
target and a length of at least one primer or probe complimentary to the target. 

• The hybridization thermodynamics may be calculated for a plurality 
of primers or probes complimentary to the target, 

10 The hybridization information may represent at least a section of a 

target and a primer or probe. 

A length of the target may be longer than a length of the primer or 
. probe and wherein the hybridization thermodynamics are calculated for a best 
target/primer or target/probe complex and for competitive mismatch complexes. 

15 Hybridization information may represent at least a section of a target 

and a primer or probe and wherein a length of a target is longer than the length of 
the primer or probe and wherein the hybridization thermodynamics are calculated 
for a best target/primer or target/probe complex and for competitive target/primer 
or target/probe complexes. 

20 The program may further execute the step of calculating 

concentration of each species in a solution at a plurality of temperatures. 

Hybridization information may also represent a primer or probe and 
wherein the length of the target is longer than a length of the primer or probe and 
wherein the hybridization thermodynamics are calculated for a best target/primer 
25 or target/probe complex and for competitive mismatch complexes and wherein the 
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fc*mtiT - ^ h * J# - program may execute the step of calculating concentration of every species in a 

solution at a plurality of temperatures. 



The hybridization thermodynamics may be calculated for at least two 
best target/primer or target/probe complexes and for their corresponding 
5 competitive mismatch complexes and wherein the program may execute the step of 
correcting for any interactions between the at least two best target/primer or 
target/probe complexes and their components. 



The above objects and other objects, features, and advantages of the 
present invention are readily apparent from the following detailed description of the 
10 best mode for carrying out the invention when taken in connection with the 
accompanying drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIGURE 1 is a schematic drawing wherein multiple equilibria are 
considered for concentration calculations; 

15 FIGURE 2a is a schematic drawing of a user input interface wherein 

the user provides various input information for a first module of the invention; 

FIGURE 2b is a schematic drawing of a user output interface wherein 
a computer provides output information corresponding to the input information of 
Figure 2a; 



20 FIGURE 3a is a schematic drawing of a user input interface wherein 

the user provides various input information for a second module of the invention; 

FIGURE 3b is a schematic drawing of a user output interface wherein 
a computer provides output information corresponding to the input information of 
Figure 3a; 
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FIGURE 4a is a schematic drawing of a user input interface wherein 
the user provides various input information for a third module of the invention; 

FIGURE 4b is a schematic drawing of a user output interface wherein 
a computer provides output information corresponding to the input information of 
5 Figure 4a; 

FIGURE 5a is a schematic drawing of a user input interface wherein 
the user provides various input information for a fifth module of the invention; 

FIGURE 5b is a schematic drawing of a user output interface wherein 
a computer provides output information corresponding to the input information of 
10 Figure 5a; 

FIGURE 6 is a block diagram flow chart illustrating the solution of 
conservation equations of the present invention; 

FIGURE 7 is a schematic diagram illustrating multiplex PCR design; 

FIGURE 8 shows prediction of molecular beacon net hybridization 
15 thermodynamics; 

FIGURE 9 shows simulation of molecular beacon hybridization 
concentrations at temperatures from 0 to 100°C; 

FIGURE 10 is a diagram of match vs. mismatch hybridization; 

FIGURE 11 shows match vs. mismatch hybridization simulation at 
20 different temperatures; 

FIGURE 12 shows a general case of competitive hybridization 
equilibria that can be solved using the described numerical methods; and 
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FIGURE 13 is an example of simultaneous equations for the general 
five molecule case. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

In general, the method and system of the present invention include 
5 rigorous thermodynamic parameterization for Watson-Crick base pairs, internal 
mismatches, terminal mismatches, terminal dangling ends, co-axial stacking 
interactions, sodium and magnesium salt dependence, denaturants (urea, formamide, 
DMSO). In addition, loop parameters for hairpins, internal loops, bulges, and 
multibranched loops are included. For DNA essentially all the parameters have 

10 been previously published or all included in the Appendix hereto. Specifically, the 
parameters which have been published include Watson-Crick parameters, sodium 
dependence, GT, GA, CT, AC, AA, CC, GG, and TT mismatches. The 
parameters included herein include dangling ends, terminal mismatches, DNA loop 
parameters, and co-axial stacking parameters. For RNA, the parameters have been 

15 published by Douglas H. Turner et al. For DNA/RNA hybrid duplexes, the 
parameters have been published by Naoki Sugimoto. 

The method and system are adapted for future implementation of 
parameters for modified nucleosides (including but not limited to inosine, 5- 
nitroindole, PNA, MOE-modified RNA, and iso-bases). With these parameters, it 

20 is possible to predict the melting temperature, Tm, of a duplex within 2°C on 
average. Correction for surface effects for DNA chip arrays is also implemented. 
In addition to predicting duplex hybridization, the software accounts for single- 
strand secondary structure. This is accomplished by a new numerical procedure for 
solving complex coupled equilibria (multi-state model). With this approach, it is 

25 possible to accurately predict not only the Tm for hybridization but also the 
concentration of every species in the solution (e.g. match duplex, mismatch, duplex, 
folded target, folded primer, primer dimer, etc.) at every temperature from 0 to 
100°C. Thus, it is possible to use this software to design oligonucleotide 
hybridization with optimized temperature, salt, and strand concentrations. 
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Predicting Accurately Primer Target Interaction Stability 

The stability of a primer/target or probe/target complex can be 
described by the free energy of association of the probe and the target. The most 
accurate way to calculate free energy of association is to use the nearest-neighbor 
5 model with accurate thermodynamic parameters. Thermodynamic parameters 
should account for Watson-Crick base pairs (SantaLucia et al., 1996; Allawi & 
SantaLucia, 1997; SantaLucia, 1998), single mismatches (Allawi & SantaLucia, 
1997), terminal mismatches (disclosed herein), dangling ends (Bommarito, Pugret 
& SantaLucia, 2000) and possibly double mismatches. Proper calculation of the 

10 monovalent and divalent salt dependence is also important (SantaLucia, 1998). 
Other loop motifs for hairpins, bulge, internal loops and multi-branched loops are 
important for single strand secondary structure prediction, but are often very 
crudely approximated. Moreover, when primer and target folding can occur, a set 
of coupled equilibria should be used to model the system. The nearest-neighbor 

15 model needs to be used to determine the equilibrium constant of each equilibrium. 
The determination of possible primer or target folding can be addressed by using 
secondary-structure prediction algorithms like M. Zuker's MFOLD (Zuker, 1989). 

Secondary Structure and Net Hybridization Thermodynamics 

Species Concentration Calculations 

20 Consider a system of strands SI and S2 with four states: folded 

target, folded probe, probe bound to target, and random coil target and probe. The 
model can be described by three equilibria as shown in Figure 1. 

The concentrations of every species for such a system can be 
analytically determined. The three equilibrium constants for such a system are 
25 shown below: 

si+s2 -" dh » 
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S2=H2 H-ra (3) 



where SI, S2, HI, H2, and DH are the random coil SI, the random coil S2, the 
folded strand HI, the folded strand H2, and the double helix DH, respectively. The 
5 conservation of SI and S2 leads to the following equations: 

C™ = SU H1+ DH (4) 



Clf = S2+H2+DH (5) 



where Csf* 1 are the total concentrations of SI and S2. [DH] and [S2] can be 
10 expressed as a function of [SI] by substituting the [HI] obtained from Equation 2, 
in Equation 6, and then substituting the [DH] obtained by Equation 6 in Equation 
1. 

[DHj^C^-lSll-KASl] (6) 



r 521 _ im r 521 _ C™-[S1]-K 2 [S1] 



15 Substitution of [H2], [DH], and [S2] from Equations 3, 6 and 7 in 

Equation 5 leads to an expression of Kj that can be rearranged as a quadratic 
equation in [SI]: 

[SI] 2 (K, + Kj K,) + [S1](K, Ga + K 3 + K 2 K 3 - K, C S1 + K 2 + 1) - (K 3 + 1) C SI = 0 

(8) 

20 This equation is simplified by making the following substitutions: 
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a = (K,+K 2 K,) (9) 
b = (K, C™ + K 3 + K 2 K 3 - K, C™ + K 2 + 1) (10) 
c = (K 3 + 1) C™ Cll) 

The physical solution of the quadratic equation (i.e. positive root) is (Press, 1999): 



mm -»Jtt (12) 



or 



™=z^hz <13) 



The second equation has better numerical stability (Press, 1999). 
[DH], [S2], [HI], and [H2] can then be calculated using Equations 1-3. 
10 Determination of Net Free Energy 

The net free energy of hybridization is calculated as follows: 



AG°3 7 „ e ,= -*rin;^, (14) 



where 



[DH] 

[S^singlestranded ] [^single stratided ] 



15 where [SI single stranded] and [S2 single stranded] are the concentrations of SI and 
S2 either in the random coil state or the hairpin states, at the temperature of the 
simulation. Using the conservation of SI and S2, Equation 14 is rewritten as 
follows: 



-16- 



WO 01/94611 



PCT/US01/18424 



AG °^ = ~ RT In ccr - [mSSs - - [AHD (16) 

Note that AG° r m has the unusual property that it depends on the total strand 
concentrations, and C^f 1 . The net free energy expresses the duplex formation 
equilibrium free energy corrected for secondary-structure formation in the single 
5 strands. 

Determination of Net Melting Temperature 

If the strands are non self-complementary two cases have to be 
considered depending on the relative strand concentrations: 

1) If SI is the limiting reagent (Clf* < Cgf), at T M : 



10 [DH]=^C™ (17) 

The concentrations of strands [SI] and [S2] are given by the following relations: 



Clf = j c «" + K ^ + W <18a) 



[S "=2(^l) (18b) 



Cl? = \ Clf + K 3 [S2] + [52] (19a) 



^ Total }_ pTotal 

15 m. "' K f; < 19 ") 
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The replacement of [SI] and [S2] in Equation 1 gives: 



scrota! J. /-i Total 



A - Lr*T"*> - rectal , K 2 + K * K * ± ^3 + 1 (20b) 

2 AT 



Using the relation AG° X = -R T In K, Equation 20 is arranged as follows: 



VP) 



i / >AG' r (2) -Ag* r (2)-AGV(3) -AGV(3) \ AG^- 

(21) 

AG° T can then be decomposed as aG° t =AH°- T AS° (assuming AC p^o) to obtain: 

1 -A/f°(2)+AJf o (0 4AS°(2>-AS <> (1) 

0=-Cjf" + C™ + e RT e R 

-Ai/»(2)-A/7 0 (3>hA// 0 (l) AS 0 (2)*AS 0 (3)-AS°(1) -A// 0 (3)+AJ/°(l) AS a (3)-A5 a (l) A//°(l) -AS'QX 

+ e ^ e R \e RT e R + e RT e R 

(22) 

10 The above equation can be solved by bisection or other numerical techniques to find 
T. This solution is the net melting temperature. 

2) If S2 is the limiting reagent (Cg* < QV\ the following relation can 
be deduced by a similar approach: 



0 = ^Cjr - C™ + ^±^±^±i (23) 

15 Again, application of the bisection method to an equation symmetric to Equation 22 
affords the net melting temperature. 
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If the strand S is self complementary, the reactions are described by 
the following equilibria: 



S + S^DH *i = ^p (24) 



S-H K 2 = ¥± (25) 



AtT^[DH]-\cl otal (26) 
The strand conservation equation is: 

c Total = [5 j + ^ + qjjgj (2?) 

Insertion of [H] and [DH] from Equations 25 and 26 in Equation 27 leads to: 



Total 

■.Total _ oak _l i\rci~ rci ^ (28) 



10 Introduction of [S] in Equation 24 gives: 



*> - w = ~ *cr ^ ° 0 • *■ * ' * ,cr * 1 (29> 



Using the relation AG° T = -R T In K, Equation 29 is rearranged as follows: 



-2AGV(2) -AG' r (2) AGV(l) 

0=e * r +2e w - C™e * r +1 (30) 
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AG° T can then be decomposed as AG° T = AH°- T AS° to obtain: 

-2A//°(2) +2A5»(2) -AH*(2) +AS°(2) -A//°(l) +AS°Q) 

0=e RT e R +2e * T e * -C T s otal e RT e R + 1 (31) 

This equation can be solved by bisection to afford the net melting temperature. 

An experimentally validated example of the accuracy of the net 
5 hybridization thermodynamics is shown in Figure 8 for molecular beacons. At the 
top of Figure 8 are the predicted thermodynamics for simple duplex formation 
assuming no competing single strand secondary structure. Using Module 1 of the 
invention, these results are similar to what would be predicted using other 
commercial software (such as oligo 6.0), though our thermodynamic database 

10 includes the dangling end effects and salt corrections are more accurate than other 
software. The middle of Figure 8 shows the single strand folding at the molecular 
beacon as output from DNA-MFOLD. The bottom table of Figure 8 shows the 
experimentally determined A6 (effective) and Tm (effective) published in Bonnet et 
al. 1999, as well as the effective Tm and A6 (effective) predicted with Module 1 

15 using the coupled equilibria calculations. Note the close agreement between 
experiments and predictions in the bottom table and the disagreement between 
experiments and the predictions using the naive simple hybridization calculation (top 
table of Figure 8). Also note the good agreement in the bottom table for the fully 
matched A-T sequence and mismatch A-A, A-C, and A-6 sequences, thus validating 

20 the mismatch parameters. 

Further, the net hybridization calculations can be extended to 
different temperatures as shown in Figure 9, to reveal how the concentrations of all 
species change with temperature. Given the extinction coefficients and fluorescence 
quantum yields, the concentration vs. temperature profiles shown in Figure 9 can 
25 be used to calculate the fluorescence vs. temperature profile (not shown), thereby 
allowing the prediction of the temperature which produces the maximum 
fluorescence signal and minimum background fluorescence signal. 
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Another manifestation of the concentration calculations is for match 
vs. mismatch discrimination (Figure 10), whereby the concentrations of all species 
at all temperatures can be calculated (Figure 11). For the particular case shown, 
optimal match vs. mismatch discrimination is predicted to occur at 0°C. The 
5 concentration calculations can be generalized for cases in which molecules can form 
many different competing unimolecular, biomolecular, and higher order complexes 
(Figure 12) using generalized equations such as shown in Figure 13 for the five 
molecule case, and solved using the algorithm in Figure 7. 

Algorithm 

10 The hybridization prediction algorithm of the present invention is 

based on a nearest-neighbor-model analysis of the sequences. The algorithm 
accounts for structural motifs including Watson-Crick base pairs (Allawi and 
SantaLucia, 1997; SantaLucia, 1998; Sugimoto et al., 1995; Xia et al., 1998), 
single internal mismatches (Allawi and SantaLucia, 1997; Allawi and SantaLucia, 

15 1998; Allawi and SantaLucia, 1998; Allawi and SantaLucia, 1998; Kierzek et al., 
1999; Peyret et al., 1999; SantaLucia, 1998), double mismatches (Allawi and 
SantaLucia, 1997) coaxial-stacking interfaces (disclosed herein) (Walter and Turner, 
1994), terminal mismatches (disclosed herein) (Freier et al., 1986) and dangling 
ends (Bommarito et al., 2000; Freier et al., 1986). Once the motifs are identified 

20 and their thermodynamic contributions are added, the sum may be corrected for salt 
effects (sodium and magnesium) and the net hybridization is calculated when 
appropriate. 

Algorithm Functions 

A first or main module of the algorithm calculates the hybridization 
25 thermodynamics (AH°, AS°, AG° 37 , T M ) of a given duplex. Net hybridization 
accounting for secondary structure in both strands is also calculated. 
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Parameterization 

Parameters are organized in three arrays. The first array contains 
internal element parameters: Watson-Crick nearest neighbors and single mismatch 
nearest neighbors. The second array contains terminal element parameters: terminal 
5 mismatches and dangling ends. A single parameter is used to account for double 
mismatches except for tandem GT mismatches, which are explicitly enumerated 
(Allawi & SantaLucia, 1997). The third array contains coaxial-stacking parameters 
(contained herein). 

For DNA sequences, the thermodynamic contribution of all Watson- 
10 Crick nearest neighbors and single internal mismatches has been systematically 
studied (Allawi and SantaLucia, 1997; Allawi and SantaLucia, 1998; Allawi and 
SantaLucia, 1998; Allawi and SantaLucia, 1998; Peyret et al., 1999). A limited 
number of sequences containing double mismatches has also been studied (Allawi 
and SantaLucia, 1997). The contributions of dangling ends (Bommarito et al., 
15 2000) have also been systematically analyzed. Salt corrections are available for 
sodium in the range 0.01 to 1 M (SantaLucia, 1998). 

For RNA sequences, the thermodynamic contribution of all Watson- 
Crick nearest neighbors has been systematically studied (Xia et al., 1998). A limited 
number of sequences containing single mismatches has also been studied (Kierzek 
20 et al., 1999). The contribution of dangling ends and terminal mismatches has also 
been systematically analyzed (Freier et al., 1986). No salt correction has been 
developed for RNA and therefore the DNA salt corrections are assumed. These 
corrections are likely to be deficient in the case of RNA. 

For DNA/RNA hybrids, the thermodynamic contribution of all 
25 Watson-Crick nearest neighbors has been systematically studied as well as a limited 
number of sequences containing single mismatches (Sugimoto et al., 1995). As no 
salt correction has been developed for DNA/RNA hybrids, the DNA corrections are 
assumed. The applicability of these corrections to DNA/RNA hybrids has not been 
tested. The parameter arrays are designed to easily accommodate implementation 
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of new parameters and salt corrections including thermodynamics parameters for 
modified bases and denaturant effects. 

Correction for Hybridization to DNA Microchips 

A linear correction of the free energy is implemented in the algorithm 
5 of the invention to correct for hybridization to DNA microchips: 

AG° 37 (microchip) =aAG° 37 (solution) +b (32) 

where a and b are user defined real coefficients. Fotin et al. (Fotin et al., 1998) 
showed that a linear relationship could be used to relate the free energies obtained 
for hybridization in solution and on microchip surfaces. However, the relation 
10 between thermodynamics measured in solution and thermodynamics measured using 
microarrays is still unclear and appears to be different depending on the 
manufacture and type of microarrays. 

User Interface: Input and Output 

Figure 2a shows the user interface input. The users enter the 
15 sequence of each strand, the hybridization conditions (hybridization temperature, 
strand concentrations, and monovalent cations and concentrations), and 
thermodynamic corrections for single strand folding. Figure 2b shows the output 
corresponding to the input in Figure 2a. 

The algorithm can be used via the Internet at: 
20 http://islLchem.wavne.edu/Hvther/hvthermlmain.html . The algorithm may be 
written in FORTRAN 77 and run on UNIX environment or other languages and 
environments. 
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Molecular Beacons ** 

The algorithm may be used to predict the thermodynamics of a set 
of literature measurements for molecular beacons (Bonnet et al., 1999). Molecular 
beacons are high specificity probes that are efficient for mutation analysis 
5 (Giensendorf et al., 1998) and multiplex detection of single nucleotide variations 
(Marras et al., 1999). The design and efficiency optimization of these beacons is 
helped by hybridization thermodynamics prediction. Bonnet et al. studied, the 
hybridization of the molecular beacon 5 CGC, TCC, CAA, AAA, AAA, AAA, 
CCG AGC G 31 to a set of four different targets including a perfect match duplex, 
10 and three different duplexes containing one mismatch. Free energy and enthalpy for 
duplex folding may be calculated using the DNA MFOLD program 
(http://mfold2.wustl.edu/-mfold/dna/fonnl.cgi). These parameters may then 
incorporated as secondary structure corrections in Figure 2a. 

The software to implement the algorithm may be written in 
15 FORTRAN, C ++ , Visual Basic, HTML, and JAVA script computer languages. 
Two graphical user interfaces may be provided: Windows application and web 
browser format. The software may run on IBM/PC, Sun, and Silicon Graphics 
platforms. 

The software may be written in several modules as described below. 

20 A. Interactive Mode: Command Line Interface in MS-DOS 

MODULE 1 (As Previously Described abovel 

Function , Module 1 predicts the hybridization thermodynamics of a given 
duplex (DNA/DNA, RNA/RNA, or DNA/RNA). 
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Input (Figure 2a) 

Input of Sequences 

1 . Only the following characters are accepted: A, a, C, c, G, g, 
T, t, U, u, /,*,+. Single blank characters and numbers will be automatically 

5 edited, but more than one carriage return is not permitted. 

2. If the duplex contains a dangling end on a strand, the 
sequence of the other strand should contain a * at the corresponding position. (This 
is very important to include for primer binding to a large target sequence). Note: 
The top strand must be entered in 5* to 3* orientation, but the bottom strand must 

10 be entered in 3* to 5* orientation. Also, a " + " must be added at the end of each 
sequence. There is a length limit of 1024 characters for sequence entries. In 
module 1, it is important to be sure that both sequences have the same length. 

Example: AAA ACCCCTGA + 
*TTTGGGGAC* + 

15 > 3. Only the bottom strand may contain coaxially stacked 

nucleotides. A should be inserted at the site of a strand nick (i.e. between the 
coaxially stacked nucleotides). This feature is useful for predicting stacked 
hybridization stability. 

Example: AAAACCCCC + 
20 TTTT/GGGG+ 

Input of Salt and Strand Concentrations 

The monovalent salt should be the sum of all monovalent cation 
concentrations in a solution in units of molarity. For example, a solution of 100 
mM KC1, 50 mM NaCl, 10 mM Na 2 P0 4 , 0.1 mM Na^DTA would account for a 
25 total of 0. 1702 M monovalent. The thermodynamic predictions are applicable over 
a salt range of 0.01 to 1 M monovalent cation. The correction applied is from 
SantaLucia (1998) Proc. Natl. Acad. Sci. 95, 1460. The sodium correction applies 

-25- 



WO 01/94611 



PCT/US01/18424 



for oligonucleotides with fewer than about 30 base pairs. For longer duplexes a 
polymer correction is required, but this is not currently implemented. 

Strand concentrations are entered in units of molarity. The program 
will accept virtually any physically relevant strand concentration. 

5 Hybridization temperature is in Celsius degrees. The limits are 0 to 

100 degrees. 

Special corrections for single-stranded secondary structure and for 
surface corrections for hybridization arrays can be input. The units for input AG° 
are kcal/mol. To determine estimates of single-strand folding energies, see Michael 
10 Zuker's RN A or DNA-MFOLD servers (see http: 
//mfold2.wsutl.edu/-mfold/diia/forml .cgi) . The current thermodynamic prediction 
software incorporates the special corrections for single-stranded secondary structure 
and for surface corrections for hybridization arrays. 

For DNA chip arrays, a linear correction can be applied. The user 
15 inputs the slope and intercept coefficients. Based on the work of Mirzabekov 
group, a slope of +1.1 and intercept of +3.2 are appropriate (see Fotin et al. 
(1998) Nucleic Acids Res. 26, 1515-1521). 

Output (Figure 2b) 

Module 1 outputs the hybridization thermodynamics at 1.0 M NaCl 
20 and 37 °C (the conditions under which the thermodynamic predictions are most 
accurate), under the salt temperature conditions specified by the user, and also 
displays the net hybridization Tm and AG 0 if the user specifies that special 
corrections are needed (this allows for single-strand secondary structure of both the 
target and probe DNA to be accounted and for surface effects of chip arrays). 
25 Predictions of AG°, AH°, AS°, and Tm are provided. 
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MODULE 2 



Function . Module 2 finds the best primers of given length 
complementary to a long target nucleic acid. DNA/DNA, RNA/RNA, DNA/RNA 
hybridization types are accepted. The user selects the number of primers to output, 
5 and the program finds the most stable primers and gives their hybridization position 
and thermodynamics of each primer. 

Invut (Fizure 3a) 



Input of Salt and Strand Concentrations 

The input of strand and salt concentrations is similar to Module 1. 

i 

10 Input of Sequences 

The target sequence is input as in Module 1 . 

Output (Figure 3b) 



Primer Length and Number of Best Primers 
Module 2 displays "number of best primers" best primers of length 
15 "primer length* in order of decreasing stability . 

Output 

Module 2 outputs "number of best primers" best primers of length 
"primer length" in order of decreasing stability along with their hybridization 
thermodynamics. 



20 MODULE 3 



Function . Module 3 walks a given primer along a given target and 
finds the thermodynamics for the best target/primer complex and for the competitive 
target/primer complexes: DNA/DNA, RNA/RNA, DNA/RNA, hybridization types 
are accepted. 
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Input (Figure 4a) 

Input of Sequences 

The input is similar to Module 1. The target has to be longer than 

the primer. 

5 Input of Salt and Strand Concentrations 

The input of salt and strand concentrations is similar to Module 1. 

Percent Stability p of Alternative Binding Sites 
Compared to the Most Stable Binding Site 

This parameter excludes all competitive sites that are not within the 

10 defined percent of the best primer stability. If the best primer stability is -5 

kcal/mol and p= 10 then any competitive site of energy higher than -5 + (10/100*5) 

= -4.5 kcal/mol will not be displayed. 

Number of Base Pairs Required to Compute the Solution 
This parameter excludes all competitive sites that contain less 
15 Watson-Crick base pairs than the defined value. 

Output (Figure 4b) 

Module 3 outputs the best primer binding site and the competitive 
binding sites that pass the filtering criteria (percent stability p of alternative binding 
sites compared to the most stable binding site and number of best primers). 

20 MODULE 4 

Function . Batch mode calculations (see below). 
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MODULE 5 

Function . Module 5 is a combination of Modules 2 and 3 and finds 
the n best primers of given length complementary to a given section of a target and 
display the thermodynamics of the target/primer system(s). Then, each best primer 
5 is walked along the whole target to find the competitive hybridization sites. The 
thermodynamics of the target/primer systems at these alternative sites is then 
displayed. DNA/DNA, RNA/RNA, DNA/RNA, hybridization types are accepted. 

Input (Figure 5a) 

Input of Sequences 
10 The target sequence is input as in Module 1. 

Input of Salt and Strand Concentrations 

The input of salt and strand concentrations is similar to Module 1 . 

Sequence Section Where to Find the Best Primers 
Module 5 finds the best primers in the target region ranking from 
15 "position of initial nucleotide" to "position of final nucleotide" . Note that Module 
5 then looks for competitive sites of each best primers in the whole target. 

Percent Stability of Alternative Binding Sites 
Compared to the Most Stable Binding Site 

The function of this parameter is the same as in Module 3. This 

20 parameter is input for each best primer corresponding to the "number of best 

primer" specified. 

Number of Base Pairs Required to Compute the Solution 
The function of this parameter is the same as in Module 3. This 
parameter is input for each best primer corresponding to the "number of best 
25 primer" specified. 
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OutpuhFieure 5b) 

Primer Length and Number of Best Primers 
Module 5 displays "number of best primers " best primers of length 
"primer length" by order of decreasing stability. 

5 Output 

Module 5 displays "number of best primers'' best primers and their 
competitive sites by order of stability along with their hybridization 
thermodynamics. The best primer and its ranked competitive hybridization sites are 
listed first. Then, the second best primer is listed with its competitive hybridization 
10 sites. 

MODULE 6 

Function . Module 6 is similar to Module 3 and walks a given primer 
along a given target and finds the thermodynamics for the best target/primer 
complex and for the competitive target/primer complexes: DNA/DNA, RNA/RNA, 
15 DNA/RNA, hybridization types are accepted. Then, Module 6 simulates the 
concentration of every species at every degree from 1 to 100 °C, as illustrated in 
Figure 6. 

Input (Not Shown) 

Input of Sequences 

20 The input is similar to Module 1. The target has to be longer than 

the primer. 

Input of Salt and Strand Concentrations 

The input of salt and strand concentrations is similar to Module 1. 
Percent Stability p of Alternative Binding Sites 
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Compared to the Most Stable Binding Site 

This parameter excludes all competitive sites that are not within the 
defined percent of the best primer stability. If the best primer stability is -5 
kcal/mol and p=10, then any competitive site of energy higher than: -5 4- 
5 (10/100*5) = -4.5 kcal/mol will not be displayed. 

Number of Base Pairs Required to Compute the Solution 
This parameter excludes all competitive sites that contain less 
Watson-Crick base pairs than the defined value. 

Correction for Target/Target Interaction, Target folding, 
10 Primer/Primer Interaction and Primer Folding 

The user is asked if he wants to correct for the interactions above. 

If the answer is "y", the user is prompted for AH° 37 corresponding to the 

interaction. Secondary structure thermodynamics can be determined using the 

Zuker algorithm as discussed in Module 1 section. 

15 Output (Not Shown) 

Concentration Output Filename 

The results from the concentration simulations (concentration of 
species at every temperature) are saved in this file. 

Output 

20 Module 6 outputs the best primer binding site and the competitive 

binding sites that pass the filtering criteria (percent stability p of alternative binding 
sites compared to the most stable binding site and number of best primers). The 
• concentration simulations are saved in a file specified by the user. 
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MODULE 7 



Function . Module 7 is a combination of Modules 2 and 5 and finds 
the n best primers of given length complementary to a given section of a target and 
display the thermodynamics of the target/primer system(s). Then, each best primer 
5 is walked along the whole target to find the competitive hybridization sites. The 
thermodynamics of the target/primer systems at these alternative sites is then 
displayed. DNA/DNA, RNA/RNA, DNA/RNA hybridization types are accepted. 
Then, Module 7, like Module 6, simulates the concentration of every species at 
every degree from 1 to 100°C, as illustrated in Figure 6. 

10 Input (Not Shown) 

Input of Sequences 

The target sequence is input as in Module 1 . 
Input of Salt and Strand Concentrations 

The input of salt and strand concentrations is similar to Module 1. 



15 Sequence Section Where to Find the Best Primers 

Module 7 finds best primers in the target region ranking from 
"position of initial nucleotide" to "position of final nucleotide. " Note that Module 
7 then looks for competitive sites of each best primers in the whole target. 



Percent Stability of Alternative Binding Sites 
20 Compared to the Most Stable Binding Site 

The function of this parameter is the same as in Module 3. This 

parameter is input for each best primer corresponding to the "number of best 

primer" specified. 



-32- 



WO 01/94611 



PCT/US01/18424 



Number of Base Pairs Required to Compute the Solution 

The function of this parameter is the same as in Module 3. This 

parameter is input for each best primer corresponding to the "number of best 

primer" specified. 

5 Correction for Target/Target Interaction, Target folding, 
Primer/Primer Interaction and Primer Folding 

For each best primer, the user is asked if he wants to correct for the 

interactions above. If the answer is *y", the user is prompted for AK° and AG° 37 

corresponding to the interaction. Secondary structure thermodynamics can be 

10 determined using the Zuker algorithm as discussed in Module 1 section. 

Concentration Output Filenames 

For each best primer, the results from the concentration simulations 
(concentration of species at every temperature) are saved in this file. The user has 
to select a different filename for each best primer. 

15 Output (Not Shown) 

Output 

Primer Length and Number of Best Primers 

Module 7 displays "number of best primers" best primers of length 
"primer length" by order of decreasing stability. 

20 Module 7 displays "number of best primers" best primers and their 

competitive sites by order of stability along with their hybridization 
thermodynamics. The best primer and its ranked competitive hybridization sites are 
listed first. Then, the second best primer is listed with its competitive hybridization 
sites. For each best primer, a file named by the user contains the concentration 

25 simulations. 
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Module 7 allows the user to design optimal primers for applications^ 
where multiple simultaneous hybridization reactions are occurring, including match 
vs. mismatch hybridization, molecular beacons, DNA oligonucleotide arrays, and 
multiplex PCR. 

5 One commercially important example for the use of Module 7 for 

primer design in a complex hybridization solution is Multiplex PCR, as shown in 
Figure 7. Module 7 allows the user to design optimal primers for Multiplex PCR 
where multiple primers have equal stabilities in binding to the target DNA. Several 
primers must be designed to specifically bind to different sites on target DNA at a 
10 given temperature with minimal background binding to mismatch sites and with 
minimal cross-hybridization between pairs of primers. 

Module 7 minimizes potential primer dimer formation and mismatch 
hybridization for all combinations of input primers. Module 7 optimizes primer 
sequence position, length, and concentration for each primer in relation to all other 
15 species in solution and provides a hybridization profile at all temperatures from 0 
to 100°C. 

Batch Mode 

MODULE 4 

Function . Module 4 allows any of the previous modules to be run 
20 in batch mode using text files to submit the input and having the data output as text 
files also. 

Type of Input Files 

There are two types of input files: 1) parameter input file, and 2) 
sequence input file. Parameter input files describe what modules to run with what 
25 hybridization parameters and on how many sequences to run them. Example of 
parameter input files for each module with comments are given in the "Batch mode 
parameter files folder. " Sequence files contain the sequences that are going to be 
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hybridized in the conditions described by the parameter input files. Examples of 
parameter input files for each module with comments are given in the "Batch mode 
sequence files folder. * 

Note that a parameter file can successively run different modules on 
5 various different sequences. 

The user is successively asked for the names of the parameter input 
file, the sequence input file and the thermodynamic data output file. Note that these 
files have to be in the directory containing the executable version of the software. 
Output files will also be created in this same directory. Names of the concentration 
10 simulation files are specified in the parameter input files. 

Examples of Batch Mode Parameter Files 

Comments in parentheses describe the meaning of each entry (note 
that an actual parameter file must not contain these comments). 

DUP (Module 1: Simple duplex calculations) 

15 1 (Number of sequences to apply this parameter file to) 

1 (Monovalent cations concentration mol/L)) 

1 Mg 2+ concentration mol/L) 

37.0 (Hybridization temperature) 

4e-4 (Top strand concentration mol/L) 

20 4e-4 (Bottom strand concentration mol/L) 

1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 

0 (Correction for top strand folding: AG° 37 ) 

0 (Correction for top strand folding: AH° 37 ) 

25 0 (Correction for bottom strand folding: AG° 37 ) 

0 (Correction for bottom strand folding: AH° 37 ) 

END (End of file required) 
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NBP (Module 2: N-best primers) 

1 (Number of sequences to apply this parameter file to) 

1 (Monovalent cations concentration mol/L) 

1 (Mg 2 * concentration mol/L) 

5 37 (Hybridization temperature) 

4e-4 (Top strand concentration mol/L) 

4e-4 (Bottom strand concentration mol/L) 

4 (Primer Length) 

3 (Number of best primers) 

10 1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 

END (End of file required) 



PWA (Module 3 primer walk match vs. mismatch sites identification) 

1 (Number of sequences to apply this parameter file to) 

15 1 . (Monovalent cations concentration mol/L) 

1 (Mg 2+ concentration mol/L) 
37 (Hybridization temperature) 
4e-4 (Top strand concentration mol/L) 
4e-4 (Bottom strand concentration mol/L) 

20 90 (Percent window of best primer stability for alternative sites) 

2 (Number of WC base pairs required to compute the solution) 
cgcg -I- (Primer sequence, + required) 

1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 

25 END (End of file required) 



BPW (Module 5 displays "number of best primers" best primers and their 
competitive sites by order of stability along with their hybridization 
thermodynamics) 

1 (Number of sequences to apply this parameter file to> 
30 1 (Lower limit of primer search area) 

10 (Upper limit of primer search area) 
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1 Mg 2+ concentration mol/L) 

37 (Hybridization temperature) 

4e-4 (Top strand concentration mol/L) 

4e-4 (Bottom strand concentration mol/L) 

5 4 (Primer length) 

1 (Number of best primers) 

1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 

800 (Percent window of best primer stability for alternative sites) 

10 2 (Number of WC base pairs required to compute the solution) 

END (End of file required) 

PWC (Module 6 primer walk with concentration calculations) 

1 (Number of sequences to apply this parameter file to) 
1 (Monovalent cations concentration mol/L) 

15 1 Mg 2+ concentration mol/L) 

37 (Hybridization temperature) 

4e-4 (Top strand concentration mol/L) 

4e-4 (Bottom strand concentration mol/L) 

90 (Percent window of best primer stability for alternative sites) 

20 2 (Number of WC base pairs required to compute the solution) 

cgcg+ (Primer sequence, + required) 

1 (Correction for microchips: slope) 

0 (Correction for microchips: intercept) 

0 (Correction for target folding: AG° 37 ) 

25 0 (Correction for target folding: AH 0 ) 

0 (Correction for target/target interaction: AG° 37 ) 

0 (Correction for target/target interaction: AH°) 

0 (Correction for primer folding: AG° 37 ) 

0 (Correction for primer folding: AH°) 

30 0 (Correction for primer/primer interaction: AG° 37 ) 

0 (Correction for primer/primer interaction: AH°) 

outconc (Concentration output file name) 
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END (End of file required) 





BWC 


(Module 7, N-best primers, primer walk, and concentration 






calculations) 




1 


(Number of sequences to apply this parameter file to) 


5 


1 


(Lower limit of primer search area) 




10 


(Upper limit of primer search area) 




1 


(Monovalent cations concentration mol/L) 




0 


Mg 2+ concentration mol/L) 




37 


(Hybridization temperature) 


10 


4e-4 


(Top strand concentration mol/L) 




4e-4 


(Bottom strand concentration mol/L) 




4 


(Primer length) 




1 


(Number of best primers) 




1 


(Correction for microchips: slope) 


15 


0 


(Correction for microchips: intercept) 




800 


(Percent window of best primer stability for alternative sites) 




2 


(Number of WC base pairs required to compute the solution) 




0 


(Correction for target folding: AG° 37 ) 




0 


(Correction for target folding: AH°) 


20 


0 


(Correction for target/target interaction: AG° 37 ) 




0 


(Correction for target/target interaction: AH°) 




0 


(Correction for primer folding: AG° 37 ) 




0 


(Correction for primer folding: AH°) 




0 


(Correction for primer/primer interaction: AG° 37 ) 


25 


0 


(Correction for primer/primer interaction: AH°) 




outconc 


(Concentration output file name) 




FND 






PPW 


(Module 8: walk a primer along itself to find interaction sites. PWA 






with probe = primer) 


30 


1 


(Number of sequences to apply this parameter file to) 




1 


(Monovalent cations concentration mol/L) 
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1 


Mg 2+ concentration mol/L) 


37 


(Hybridization temperature) 


4e-4 


rPHrripr rnnrf*ntT;itifvn mnl/T ^ 


900 


(Percent window of best primer stability for alternative sites) 


2 


(Number of WC base pairs required to compute the solution) 


END 


(End of file required) 



Examples of Batch Mode Sequence Files 

For Module 1: Pup 
1 (Sequence number) 

10 agcgca+ (Top strand sequence) 
tcgcgt+ (Bottom strand sequence) 

For Module 2: NBP 

1 (Sequence number) 

agcgca+ (Target sequence) 

15 For Module 3: 

1 (Sequence number) 

cgcctgcggccc+ (Target sequence) 

For Module 5: bpw 
1 (Sequence number) 

20 cgcctgcgccc+ (Target sequence) 

For Module 6: pwc 

1 (Sequence number) 

agcgca+ (Target sequence) 

For Module 7: bwc 
25 1 (Sequence number) 

agcgca+ (Target sequence) 
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For Module 8: ppw 
. 1 (Sequence number) 

agcgca+ (Primer sequence) 

Example of Batch Mode Parameter and Sequence Files 
5 to Run Different Modules Successively 



Parameter File : 

DUP (Executes Module 1) 

2 (Apply to Module 1 to 2 sequence sets) 

0.05 
10 1.5e-3 
37.0 
le-6 
2e-7 
1 

15 0. 

-2.12 
-37.3 
0 
0 

20 PWC (Executes Module 6) 

1 (Apply to Module 6 to 1 sequence set) 

0.16 

0.0025 

37 

25 10e-9 
le-9 
800 
8 

TCGAACGTAC + 
30 1 
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0 
0 
0 
0 

5 0 
0 
0 
0 
0 

10 outwash 

DUP (Executes Module 1) 

4 (Apply to Module 1 to 4 sequence sets) 

1 

0 

15 37.0 
le-6 
le-6 
1 
0 

20 0 
0 
0 
0 

END 

25 Other modules can be similarly appended. 
Sequence File 

1 (input for Module 1) 

ttgcctaggggaccaggtccaact+ 
aacggatcccctggtccaggttga+ 
30 2 (input for Module 1) 
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ttgcctaggggaccaggtccaact + 
aacggatcccctggtccaggttga + 
3 

CAGCTTGC ATGA AA AGCTTGCGTGT + 
5 4 

AAAAAA+ 

TTTTTT+ 

5 

acgcgc+ 
10 tgcgcg+ 
6 

gggaaagggg+ 

*cctttccc*+ 

7 

15 tttaaattt+ 
aaatttaaa+ 
8 

cgcgtgagggcc+ 
gcgctctccccgg+ 



(input for Module 6) 
(input for Module 1) 

(input for Module 1) 

(input for Module 1) 

(input for Module 1) 

(input for Module 1) 



20 Parameterization of the Algorithm of the Invention 



Caution: RNA/RNA and DNA/RNA duplexes, contain motifs for 
which no literature data are available. In these cases, DNA/DNA parameters are 
assumed. Therefore, predictions might be inaccurate. Users are encouraged to use 
this program with caution and discernment. 



25 No data are available for the following motifs: 

RNA/RNA single mismatches 
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RNA/DNA single mismatches 
dangling ends 
terminal mismatches 
Single mismatches 

5 Double mismatch parameters are estimated for all types of duplexes 

(DNA/DNA, RNA/RNA, DNA/RNA). 





DNA THERMODYNAMIC PARAMETERS 




Watson-Crick nearest-neiehbors 






10 


SantaLucia, Allawi, and Seneviratne 
(1996) Biochemistry 35, 3555; 

/VlldWl ana odllULUUCla \LzJy 1 ) 

Biochemistry 36, 10581 


12 parameters 


108 sequences 




Single mismatch nearest-neiehbors 






15 


AUoWl oilU. 0£UlUXJL«UUia ^ 1 77 / J 

Biochemistry 36, 100581; 
Allawi and SantaLucia (1997) 
Nucleic Acids Res., 26, 2694; Peyret 
et al., (1999) Biochemistry 38, 3468 


44 parameters 180 sequences 

Allawi and SantaLucia (1998) Biochemistry 37, 2170 
Allawi and SantaLucia (1998> Biochemistry 37, 9435 


20 


Terminal Mismatch nearest- 
neighbors 

(Appendix) 


48 parameters 


48 sequences 




Dangling end nearest-neiehbors 






25 


S. Bommarito, Peyret, SantaLucia 
(2000) Nucleic Acids Res. 28, 1929- 
1934. 


16 parameters 


16 sequences 




Na + dependence law 








SantaLucia (1998) Proc. Natl Acad. 
Sci. USA 95, 1460-1465 


1 parameter 


86 sequences 
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RNA THERMODYNAMIC PARAMETERS 


Watson-Crick nearest-neiehbors 


12 parameters 


Xia et aL (1998) Biochemistry 37, 14719 


Single Mismatches 


44 parameters 


Kierzek, Burkard, and Turner (1999) 
Biochemistry 3o, 14zl4-l4z23 


Terminal mismatch nearest-neiehbors 


48 parameters 


Freier et al. (1986) Proc. Natl Acad. ScL USA 83, 
9373 


Dangling end nearest-neiehbors 


32 parameters 


Freier et al. (1986) Proc. Natl. Acad. ScL USA 83, 
9373 


Coaxial stacking nearest-neighbors 


16 parameters 


Walter and Turner (1994) Biochemistry 33, 12715 


Loop parameters 

Matthews et al. (1999) J. Mol Biol. 288, 911-940 





HYBRID DNA/RNA THERMODYNAMIC PARAMETERS 


Watson-Crick nearest-neighbors 

Sugimoto et al. (1995) Biochemistry 34, 11211 


17 parameters 


rU^dG and rG-dT mismatches 




Sugimoto et al. (1997) Nucleic Acids Symp. Or. '37, 199 
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DNA LOOP THERMODYNAMIC PARAMETERS 



Hairpins 

Hilbers et al. (1985) Biochie 67, 685-695 

Blommers et al. (1989) Biochemistry 28, 7491-7498 

Antao et al. (1991) Nucleic Acids Res. 19, 5901-5905 

Antao et al. (1991) Nucleic Acids Res. 20, 819-824 

Senior et al. (1988) Proc. Natl Acad. Sci. USA 85, 6242-6246 



10 



Bulges 

LeBlanc and Morden (1991) Biochemistry 30, 4042-4047 
Zieba et al. (1991) Biochemistry 30, 8018-8026 
Ke et al. (1995) Biochemistry 34, 4593-4600 
Turner, D.H. (1992) Can. Opin. Struc. Biol. 2, 334-337 



15 



Multibranched Loops 

Kadrmas et al. (1995) Nucleic Acids Res. 23, 2122 
Lilley and Hallam (1984) 7. Mol. Biol. 180, 179-200 
Lu et al. (1991) J. Mol Biol. 223, 781-789 
Ladbury et al. (1994) Biochemistry 33, 6828-6833 
Leontis et al. (1991) Nucleic Acids Res. 19, 759-766 



20 



25 



The parameters for multibranched loops are from a best fit analysis of 
secondary structure predictions vs. experiments as done by Jaeger et al. for 
RNA (Jaeger et al. (1989) PNAS 86, 7706-7710). The current parameters for 
multibranched loops neglect the sequence and complicated length dependence 
described by Leontis and coworkers, but approximate 4-way junctions fairly 
well. Implementation of more complicated rules will require modification of 
the MFOLD algorithm. 



While the best mode for carrying out the invention has been 
described in detail, those familiar with the art to which this invention relates will 
recognize various alternative designs and embodiments for practicing the invention 
as defined by the following claims. 
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APPENDIX*" 
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Table 1: Thermodynamic Parameters for Duplex Formation in 1M NaCl a . 



ATGAGCTCAA 
^ACTCGAGT A 


AH ob 
(kcal / mol) 

-56.7 ± 2.5 


AS ob 
(cal/molK) 

-155.8 ± 3.7 


agV 

(kcal / mol) 
-9.07 ±0.12 


(°C) 
52.3 


AAGAGC TCTA 
ATCTCGAGAA 


-55.1 ± 1.4 


-149.2 ± 3.8 


-8.91 ± 0.12 


56.0 


AGTAGCTACA 
ACATCAATGA 


-60.3 ± 1.8 


-165.8 ± 5.1 


-9.03 ±0.12 


54.5 


ACGATATCGA 
AGCTAT A G C A 


-67.8 ± 1.3 


-192.0 ± 3.5 


-8.87 ± 0.06 


49.4 


£ T GAGC T C A £ 
£ A C T C G A G T £ 


-50.6 ± 1.2 


-136.9 ± 2.8 


-8.15 ± 0.11 


52.9 


£A GAGC T CT £ 
£TCTCGAGA£ 


-51.4 ± 1.3 


-137.5 ± 3.2 


-8.34 ±0.14 


56.9 


£ G T A G C T AC £ 
£C AT C G AT G£ 


-55.8 ± 1.8 


-154.6 ± 4.8 


-8.04 ±0.16 


49.7 


£CGATATCG£ 
£ G C T A T AGC£ 


-59.8 ± 1.2 


-166.9 ± 3.0 


-8.12 ± 0.07 


49.7 


GTGAGCTCAS 
GACTCGAGTfi 


-52.6 ± 1.3 


-140.7 ± 3.3 


-8.57 ± 0.13 


57.7 


£ A GAGC TCT£ 
£TCTCGAGA£ 


-52.3 ± 1.4 


-142.5 ± 4.1 


-8.34 ± 0.08 


52.3 


£G T AGC T AC£ 
£C ATCGATG£ 


-59.2 ± 1.0 


-164.2 ± 2.8 


-8.68 ± 0.07 


51.2 


£CGAT ATCGfi 
£GCT AT AGC£ 


-65.7 ± 0.9 


-183.6 ± 2.2 


-8.81 ± 0.07 


52.2 
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Table 1: Continued 3 . 







AH ob 
(kcal / mol) 


AS ob 
(cal/molK) 


AG° 3 - b 
(kcal / mol) 


(°C) 


I T G A G C 
I A C T C G 


T C A I 
A G T X 


-55.4 ± 1.1 


-149.4 ± 3.0 


-8.62 ±0.10 


56.9 


I A G A G C 
I T C T C G 


T C T X 
A G A X 


-56.5 ± 1.3 


-154.3 ± 3.7 


-8.72 ± 0.09 


54.3 


X G T A G C 
1 C A T C G 


T A C X 
A T G X 


-63.8 ± 0.8 


-178.0 ± 2.1 


-8.75 ± 0.06 


52.1 


X C G A T A 
X G C T A T 


T C G X 
A G C X 


-66.8 ± 0.6 


-187.9 ± 1.6 


-8.60 ± 0.04 


50.5 


£ T G A G C 
A A C T C G 


T C A A 
A G T £ 


-53.6 ± 1.3 


-145.8 ± 4.0 


-8.42 ± 0.06 


53.6 


AT G A G C 
CACTCG 


T C A £ 
A G T A 


-54.0 ± 1.3 


-144.4 ± 3.2 


-8.92 ± 0.15 


58.8 


CGTAGC 
A C A T C G 


T A C A 
AT G £ 


-56.8 ± 1.4 


-155.6 ± 3.5 


-8.53 ± 0.13 


53.7 


A G T A G C 
£ C A T C G 


T A C £ 
AT G A 


-57.1 ± 1.2 


-156.0 ± 3.0 


-8.71 ± 0.16 


54.4 


£ C G A T A 
A G C T A T 


T C G A 
A G C £ 


-61.8 ± 0.5 


-172.7 ± 1.4 


-8.30 ± 0.03 


50.6 


A C G A T A 
£ G C T A T 


T C G £ 
A G C A 


-58.3 ± 1.7 


-158.5 ± 4.2 


-8.91 ± 0.17 


56.8 


£ A G A G C 
AT C T C G 


T C T A 
A G A £ 


-54.6 ± 0.6 


-147.9 ± 1.4 


-8.66 ± 0.07 


55.4 
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Table 1: Continued 2 . 













AH ob 


AS ob 


AG° 37 b 


T M C 












(Kcai / moi^ 


{Cai / II1U1 XV/ 




(°C) 


A A 


G A G C 


T 


C T 


c 


-55.5 ± 1.4 


-150.4 ± 4.4 


-8.85 ± 0.08 


55.8 


C T 


C T C G 


A 


G A 












T T 


G A G C 


T 


C A 


£ 


-52.2 ± 1.0 


-140.8 ± 2.4 


-8.38 ± 0.11 


54.9 


£ A 


CTCG 


A 


G T 


I 












U A U t 


T 


C A 


T 


-JJ.l 3C U.o 


1<rt c + la 
-IJU.j 3: i.y 


-R 4? + ft ftR 


jj.- 


1 A 


CTCG 


A 


G T 


£ 










T G 


T A G C 


T 


A C 


£ 


-58.0 ± 1.4 


-159.6 ± 3.7 


-8.39 ±0.10 


52.9 


£ c 


A T C G 


AT G 


1 










£G 


T A G C 


T 


A C 


I 


-59.4 ± 0.9 


-164.6 ± 2.1 


-8.21 ± 0.06 


51.7 


I c 


A T C G 


A 


T G 


£ 










T C 


G A T A 


T 


C G 


£ 


-61.7 ± 1.1 


-170.6 ± 3.2 


-8.33 ± 0.12 


53.7 


£ G 


C T A T 


A 


G C 


1 










£C 


G A T A 


T 


C G 


X 


-57.9 ± 1.1 


-159.5 ± 2.7 


-8.11 ± 0.12 


52.7 






A 


G C 












1 A 


GAGC 


T 


C T 




-55.0 ± 1.4 


-148.0 ± 3.5 


-8.80 ± 0.12 


57.5 


c t 


r t r o 


A 


G A 


T 










£ A 


GAGC 


T 


C T 


l 


-51.5 ± 1.1 


-137.9 ± 2.6 


-8.44 ± 0.15 


56.4 


T T 


P T p O 


A 


G A 












G T 


GAGC 


T 


C A 


A 


-54.2 ± 1.4 


-146.1 ± 3.3 


-8.77 ± 0.16 


56.8 


A A 


C G C G 


A 


G T 


£ 










AT 


GAGC 


T 


C A 


£ 


-55.4 ± 1.2 


-148.6 ± 3.1 


-9.03 ± 0.14 


59.1 


G A 


CTCG 


A 


G T 


A 
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Table 1 : Continued 8 . 













AH ob 


AS ob 


AG° 37 b 


T c 
»M 












flrral / mol i 


(cal / mol K) 


(kcal / mol) 


(°C) 


\j Var 


TAG 

X A VJ 


c 

V-> 


T A C 


A 


-59.4 ±1.6 


-163.3 ± 3.9 


-8.76 ±0.15 


53.S 


AC 


AT C 


G 


AT G 


G 










A Car 




p 

Vv 


TAT 


VJ 


-63 7 ± 16 


-174.9 ± 4.1 


-9.46 ± 0.18 


56.4 


GC 


A T C 


G 


AT G 


A 










fi c 


n a T 
u A 1 


A 


T C CI 


* 


-60 4 4: 08 
"Ov.H VI. O 


-166 6 ± 2 1 


-8.50 ± 0.09 


53.5 


AG 


C T A 


T 


A GC 


G 












VJ xv X 


A 


T C G 

X Vv VJ 


G 


-61 1 ± 1 3 


-167.5 ± 3.4 


-9.04 ± 0.14 


55.8 


G G 


C T A 


T 


A G C 


A 










g a 

VJl <rv 


GAG 

VJ -fx VJ 


r 


T C T 

X V* x 


A 


-54.0 ±1.1 


-144.9 ± 2.7 


-8.82 ± 0.14 


57.9 


AT 


C T C 


G 


A G A 


Q 












GAG 

VJ i» VJ 


r 


T C T 

X Vy JL 


G 


-54.8 ± 1.7 


-148.0 ± 5.0 


-8.90 ± 0.09 


56.3 


<£ T 


C T C 


G 


A G A 


A 










T T 


GAG 

W VJ 


c 

Vv 


T C A 

X v<» -r» 


G 


-56 8 ± 0 7 


-155.5 ± 1.7 


-8.64 ± 0.06 


53.5 




C T C 


G 


AGT 


X 












G A H 

VJ A VJ 


c 


t r a 

X V» -TV 


T 




-156 7 ± 2 0 


-8 80 ± 0 05 


54.7 


I A 


C T C 


G 


AGT 


a 










T G 


TAG 

1 /v VJ 


r 

Vv 


T A C 

X /k V* 




_<o 2 ± l 3 


-161 3 ± 3 6 


-8 93 ± 0.08 


56.2 


QC 


AT C 


G 


AT G 


X 










G G 


TAG 


r 

V-/ 


T A C 


T 


-64 8 ± 2 6 


-182 3 ± 7 2 


-8.63 ±0.17 


50.0 


1 c 


AT C 


G 


AT G 


£ 










T C 


GAT 

VJ X 


A 


T C G 

X V** VJ 




-63 3 ± 0 6 


-176 9 ± 1 6 

— A /Vj.J F J- X .V 


-8 42 ± 0 05 


51.2 


£C 


C T A 


T 


A G C 


i 












GAT 


A 


T C G 


i 


-63.7 ± 0.9 


-177.4 ± 2.2 


-8.73 ± 0.12 


52.6 


I G 


C T A 


T 


A G C 


a 
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Table I: Continued 2 . 





AH ob 


AS 


AO 37 


T c 

* M 




(kcal / mol) 


(cal/mol K) 


(kcal / mol) 


(°C) 


X A u A u t 1 *~ 1 vj 


-J 1 .0 3: u.o 


-1 57 A ± 1 7 


-8 95 ± 0 03 


55.6 


£T CT CG AGA 1 










GAGAGCTCT T 


-57.3 ± 1.2 


-155.5 ± 3.3 


-8.95 ±0.10 


56.4 


ITCTCGAGAfi 










Core seauences 










CGATATCG d 


-51.9 ± 0.6 


-145.3 ± 1.4 


-6.89 ± 0.09 


44.1 


GCTATAGC 










GTAGC TAC d 


-51.6 ± 0.6 


-143.7 ± 1.3 


-7.01 ± 0.08 


45.6 


CATCGATG 










AGAGCTCT 


-50.0 ± 0.7 


-136.5 ± 1.7 


-7.76 ± 0.06 


50.2 


T CTCGAGA 










TGAGCTCA 


-50.5 ± 0.5 


-137.7 ± 1.3 


-7.73 ± 0.04 


50.4 


ACTCGAGT 











a The top strand of each duplex is represented in the 5' to 3' orientation and the bottom strand 
is shown in the 3 1 to 5 1 direction. Terminal mismatch nearest neighbors are represented in bold. 

Mismatches are underlined. b AH 0 , AS 0 , and AG° 37 are the error-weighted averages of the 1/T M vs. 
In C T plot and curve fit methods in Table SI. Errors reflect the precision of the data (see text). 

c T M calculated using 10" 4 M total strand concentration. d Data from reference (19). 
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Table 2: Nearest-neighbor thermodynamic paramters of 
like-with-like base terminal mismatches in 1 M NaCl 



Dimer 


AH 0 " 


AS°° 


AG° 37 


Sequence 8 


(kcal/mol) 


(e.u) 


f\ _ ^ i / » 1 1 

(kcal/mol) 




Terminal A»A Mismatches 




M/TA 


-3 1 ± 1 3 


-7.8 ± 2.0 


-0.67 ± 0.06 


TAMA 


-2.5 ± 0.8 


-6.3 ± 2.1 


-0.58 ± 0.07 


CA/GA 


-4.3 ± 1.0 


-10.7 ± 2.6 


-1.01 ± 0.07 


GA/CA 


-8.0 ± 0.7 


-22.5 ± 1.9 


-0.99 ± 0.06 




Terminal O 


>C Mismatches 




AQ/TQ 


-0.1 ± 0.6 


0.5 ± 1.5 


-0.21 ± 0.06 


TQ/AQ 


-0.7 ± 0.7 


-1.3 ± 1.8 


-0.29 ± 0.07 


CC/GC 


-2.1 ± 0.9 


-5.1 ±2.5 


-0.52 ± 0.09 


GQ/CQ 


-3.9 ± 0.7 


-10.6 ± 1.7 


-0.62 ± 0.06 




Terminal G-G Mismatches 




AS/Tfi 


-1.1 ± 0.7 


-2.1 ± 1.8 


-0.42 ± 0.07 


TQ/AQ 


-1.1 ± 0.8 


-2.7 ± 2.2 


-0.29 ± 0.05 


ca/Ge 


-3.8 ±0.6 


-9.5 ± 1.5 


-0.83 ± 0.05 


Gfi/Cfi 


-0.7 ±0.5 


-19.2 ± 1.3 


-0.96 ± 0.06 




Terminal T«T Mismatches 




AI/TI 


-2.4 ± 0.6 


-6.5 ± 1.6 


-0.45 ± 0.05 


TI/AI 


-3.2 ± 0.7 


-8.9 ± 2.1 


-0.48 ± 0.05 


CI/GI 


-6.1 ± 0.5 


-16.9 ± 1.2 


-0.87 ± 0.05 


GJ/CI 


-7.4 ± 0.4 


-21.2 ± 1.1 


-0.86 ± 0.05 



a Thermodynamic parameters and their corresponding errors are 
calculated from Table 1 using equations 4 and 5. 
b Dimers are given in antiparallel orientation (e.g. A£/TA equals 
S-ACO 1 paired with 3-TA-5'). Mismatches are underlined. 
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Table 3: Nearest-neighbor thermodynamic fftflairflfcter^jq^ 
mixed-base terminal mismatches in 1 M NaCl 



Dimer 


AH 00 


AS 0 " 


AG° 37 b 


sequence 3 


(kcal/mol) 


(e.u) 


(kcal/mol) 




Terminal A»C Mismatches 




A£/T£ 


-1.6 ± 0.7 


-4.0 ± 2.1 


-0.35 ± 0.04 


A£/TA 


-1.8 ±0.7 


-3.8 ± 1.7 


-0.59 ± 0.08 


CA/G£ 


-2.6 ± 0.8 


-5.9 ± 1.8 


-0.76 ± 0.07 


CC/GA 


-2.7 ± 0.7 


-6.0 ±1.6 


-0.85 ± 0.09 


GA/CQ 


-5.0 ± 0.4 


-13.8 ± 1.0 


-0.71 ± 0.05 


GC/CA 


-3.2 ± 0.9 


-7.1 ± 2.2 


-1.01 ±0.10 


TA/AC 


-2.3 ± 0.5 


-5.9 ± 1.1 


-0.45 ± 0.05 


TC/AA 


-2.7 ± 0.8 


-7.0 ± 2.4 


-0.55 ± 0.05 




Terminal OT Mismatches 




A£/TI 


-0.9 ± 0.5 


-1.7 ± 1.4 


-0.33 ± 0.06 


Al/T£ 


-2.3 ± 0.5 


-6.3 ± 1.2 


-0.35 ± 0.05 


C0G1 


-3.2 ± 0.8 


-8.0 ± 2.0 


-0.69 ± 0.07 


CX/GQ 


-3.9 ± 0.6 


-10.6 ± 1.2 


-0.60 ± 0.05 


G£/CI 


-4.9 ± 0.6 


-13.5 ± 1.7 


-0.72 ± 0.08 


GI/C£ 


-3.0 ± 0.6 


-7.8 ± 1.5 


-0.61 ± 0.08 


T£/AI 


-2.5 ± 0.8 


-6.3 ± 2.0 


-0.52 ± 0.07 


T£/AQ 


-0.7 ± 0.6 


-1.2 ± 1.6 


-0.34 ± 0.08 




Terminal G»A Mismatches 




AAfTQ 


-1.9 ± 0.7 


-4.4 ± 1.8 


-0.52 ± 0.08 


Afi/TA 


-2.5 ± 0.7 


-5.9 ± 1.7 


-0.65 ± 0.07 


CAGQ 


-3.9 ± 0.8 


-9.6 ± 2.1 


-0.88 ± 0.09 


V^VJ/ VXTV 


-6.0 ± 0.9 


-15.5 ± 2.1 


-1. 15 ±U.l 


GA/CG 


-4.3 ± 0.5 


-11.1 ± 1.3 


-0.80 ± 0.06 


Gfi/CA 


-4.6 ± 0.7 


-11.4 ± 1.8 


-1.08 ± 0.09 


TA7Afi 


-2.0 ± 0.7 


-4.7 ± 1.6 


-0.53 ± 0.07 


Tfi/AA 


-2.4 ± 0.9 


-5.8 ± 2.7 


-0.57 ± 0.05 
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Table 3: Continued 



Dimer 


AH° b 


AS° b 


AG° : , 7 


sequence 2 


(kcal/mol) 


(e.u) 


(kcal/mol) 




Terminal G»T Mismatches 




AG/IT 


-3.2 ± 0.4 


-8.7 ± 1.1 


-0.45 ± 0.04 


AI/Tfi 


-3.5 ± 0.4 


-9.4 ± 1.2 


-0.54 ± 0.03 


Cfi/GI 


-3.8 ± 0.7 


-9.0 ± 1.9 


-0.96 ± 0.06 


d/G2 


-66 ± 1.3 


-18.7 ±3.6 


-0.81 ± 0.09 


GS/CI 


-5.7 ± 0.4 


-15.9 ± 1.0 


-0.76 ± 0.05 


Giyce 


-5.9 ± 0.5 


-16.1 ± 1.3 


-0.92 ± 0.07 




-3.9 ± 0.5 


-10.5 ± 1.2 


-0.59 ± 0.03 


H/A2 


-3.6 ± 0.7 


-9.8 ± 1.9 


-0.59 ± 0.06 



3 Thermodynamic parameters and their corresponding errors 
are calculated from Table 1 using equations 4 and 5. 

b Dimers are given in antiparallel orientation (e.g. AQ/T& equals 
5 ! -AC-3 f paired with 3-TA-5 1 ). Mismatches are underlined. 
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Table S 1 : Thermodynamic Parameters for Duplex Formation in 1M NaCl a . 







AH 0 
(kcal / mol) 


AS 0 
(cal/mol K) 


AG 0 .,, 
(kcal / mol) 


T M b 
<°C) 


AT G A G C 

A A T" 


T C A A 

A KJ 1 & 


c -56.0 ± 


3.2 
4 0 


-151.4 ± 10.1 
-156.5 ± 4.0 


-9.07 ± 0.13 
-9.11 ± 0.46 


57.0 
56.7 


AAGAGC 

A T P T p 


TCT A 

A n A A 
A U A /\ 


c -54.6 ± 
d -5fi 0 ± 


1.8 
2 3 


-146.9 ± 
-151j6 ± 


5.3 
5.4 


-8.91 ± 0.12 
-8.97 ± 0.64 


57.7 
56.4 


A G T A G C 

A P ATP A 


T A C A 

A T P A 
A 1 lj A 


c -59.9 ± 
d 61 l ± 

-Ol.J 


2.2 


-164.0 ± 
-168.3 ± 


6.6 
7.9 


-9.03 ±0.12 
-9 09 ± 0 83 


55,4 
.55.3 


ACGAT A 

A P P T AT 1 


T C G A 

A P A 
A U L ^ 


c -60.5 ± 

d 7A 51 i 


2.4 

1 c 
i. J 


-166.5 ± 
-198.4 ± 


7.7 
3.9 


-8.86 ± 0.06 


54.3 
53 6 


£ T GAGC 
p a p t p n 


TCAC 
a n t p 

A O I ^ 


c -52.3 ± 

-JU.H 


4.0 


-142.4 ± 12.7 
-136.6 ± 2.8 


-8.16 ± 0.11 

c n< ± n 7< 
-o.uo UJO 


52.4 

S9 4 


£ A G A G C 
P T P T P (1 


T C T £ 
a n a p 

A u A t 


c -55.6 ± 

d -40 Q ± 


2.5 
i .j 


-152.1 ± 
-134.5 ± 


7.7 
3.5 


-8.36 ± 0.14 
-R 14 ± ft 4S 

-O.J*t V/.t J 


52.8 


CGTAGC 

Vv V* r\. X Vj 


T AC £ 

A T CZ P 


c -55.1 ± 
d R ± 

JU.o 


2.4 
9 7 


-151.8 * 
-156.9 ± 


7.2' 
6.6 


-8.04 ± 0.16 

-R OR db ft 7ft 


50.9 

Jtr. / 


£ C G A T A 
P n P T AT 


T C G £ 
a n p p 

A VJ L L 


c -57.3 ± 
-60.5 * 


2.6 
1.3 


-158.4 ± 
-168.3 ± 


8.1 
3.3 


-8.12 ± 0.07 

R 91 ± ft ^1 
-O.Z1 - 1 * U.j 1 


50.8 


STGAGC 


T C A £ 
A G T £ 


c -55.2 * 
4 -50.9 ± 


2.1 
1.7 


-150.4 ± 
-137.1 ± 


6.4 
3.9 


-8.58 ± 0.14 

O ^ O -L. f\ A O 

-8.38 * 0.48 


54.2 

C A A 

54.4 


£ A G A G C 
STCTCG 


T C T £ 
A G A £ 


c -56.8 ± 
d -48.0 ± 


2.0 
1.9 


-155.1 ± 
-128.1 ± 


5.6 
5.9 


-8.72 ± 0.26 
-8.30 ± O.09 


54.6 
54.9 
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Table SI: Continued. 















AH° 




AS 0 




AG 0 


37 
















(kcal / mol) 


(cal/molK) 


(kcal / mol) 


(°C) 


£G 


T A G C 


T 


A C 


£ 


C 


-56.9 




1.5 


-155.4 




4.8 


-8.67 


± 


0.07 


54.2 


G C 


A T C G 


A 


T G 


£ 




-61.1 




1.4 


-168.5 


± 


3.4 


-8.83 


± 


0.35 




£C 


G A T A 


T 


C G 


£ 


c 


-62.4 


± 


2.4 


-172.9 


± 


7.4 


-8.79 


± 


0.08 


53.3 


£G 


C T A T 


A 


G C 


£ 


d 


-66.2 




0.9 


-184.7 


± 


2.3 


-8.93 


± 


0.21 


53.0 


T T 


G A G C 


T 


C A 


I 


c 


-56.8 


± 


1.5 


-155.1 


db 


4.7 


-8.63 


± 


0.10 


54.2 


I A 


CT C G 


A 


G T 


1 


d 


7 




1 7 


-145 6 




3 9 


-8 52 


± 


0 46 


54.4 


T A 
-L A 


\J l\ VJ \-> 


T 


C T 

^ A 


T 

A. 


c 


-58.7 


± 


1.6 


-160.3 


± 


4.6 


-8.96 


± 


0.19 


55.4 


T T 


C T C G 


A 


G A 


I 


d 


-52.9 


± 


2.1 


1 A*\ C 

-142.5 


i 

db 


C A 

6.4 




± 


U.1U 


55.6 




i A VJ ^ 


T 
1 


A C 
A 


T 
A 


c 


-63.2 


± 


1.1 


-175.5 


± 


3.4 


-8.75 




0.06 




1 c 


A T C G 


A 


T G 


I 


d 


£. A A 

-64.4 




1.1 


-1 (yA 


± 


2.6 




± 


U.-iD 


52.8 


x c 


G A T A 


Hp 

T 


C G 


X 


c 


-64.6 


± 


1.3 


-180.4 


i 


4.0 


-8.59 


± 


0.05 


CI 7 

51.7 


T G 


C T A T 


A 


GC 


1 


d 


-67.4 


± 


0.7 


-189.4 


± 


1.8 


-8.69 


± 


0.16 


51.5 


r t 


G A G c 


T 


P A 


A 


c 


-56.0 


± 


2.1 


-153.1 


± 


6.2 


-8.58 


± 


0.15 




A A 


C T C G 


A 


GT 


£ 


d 


^? o 




1 7 
i. / 




-L 
-3- 




-o.->y 


4- 


n 07 

v.U / 


54.1 


<A 1 


VJ A. vJ \* 


T 




it 


c 


-57.0 


± 


2.4 


-154.9 


± 


7.1 


-8.94 


± 


0.16 




£A 


C T C G 


A 


G T 


A 


d 


*1 7 
OZ. / 




1.0 


-141. / 




J.O 


-O. / J 


•4- 




56.1 


£ G 


T A G C 


T 


A C 


A 


c 


-57.8 


± 


2.9 


-158.8 


± 


8.9 


-8.54 


± 


0.13 




AC 


A T C G 


A 


T G 


£ 


d 


-56.5 




1.6 


-155.0 




3.8 


-8.46 


± 


0.42 


53.0 


AG 


T A G C 


T 


A C 


£ 


c 


-59.4 




4.0 


-163.2 




12.2 


-8.74 


± 


0.18 


53.9 


£C 


A T C G 


A 


T G 


A 


d 


-56.8 




1.3 


-155.5 


± 


3.1 


-8.60 




0.35 


53.8 


£ C 


G A T A 


T 


C G 


A 


c 


-61.6 


db 


1.1 


-171.8 


± 


3.5 


-8.30 




0.03 


50.8 


AG 


C T A T 


A 


G C 


£ 


d 


-61.9 


± 


0.6 


-172.9 


± 


1.5 


-8.30 




0.14 


50.7 
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Table SI: Continued. 













AH° 




AS' 


) 


AG 0 


37 


»-p b 
T M 












(kcal / mol) 


(cal/mol K) 


(kcal / mol) 


(°C) 


A C 


CI A T A 

(j A 1 A 






c 


-61.8 


dt 


3.2 


-170.3 


i 


9.6 


-8.94 


X 


0.19 


54.3 


C G 


C T A T 


A G C 


A 


d 


-57.0 


± 


1.9 


-155.8 


± 


4.6 


-8.71 


4- 


0.52 


54.4 


C A 


GAGC 

VJ JM. VJ w 


T C T 


A 


c 


-55.7 




1.1 


-151.6 




3.5 


-8.67 


X 


0.07 


54.6 


AT 


C T C G 


AG A 


£ 


d 


-54.2 




0.7 


-147.2 




1.6 


-8.59 


X 


0.18 


54.6 


A A 


GAGC 


T C T 


£ 


c 


-58.6 


± 


4.1 


-159.8 


± 


12.4 


-9.07 


X 


0.26 


56.1 


£ T 


C T C G 


A G A 


A 


d 


-55.1 


± 


1.5 


-149.1 




4.7 


-8.83 


X 


0.08 


55.8 


1 T 


GAGC 


T C A 


£ 


c 


-54.1 


± 


1.9 


-147.4 


± 


5.8 


-8.40 


X 


0.12 


53.4 


£ A 


C T C G 


A GT 


I 


d 


-51.5 


± 


1.1 


-139.4 


± 


2.6 


-8.27 


X 


0.32 


53.5 


£ T 


GAGC 


T C A 


I 


c 


-57.3 


i 

X 


4.5 


-157.5 


i 

x 


14.3 


-8.43 


4- 
X 


a r\f\ 

o.oy 


52.7 


I A 


CTCG 


A G t 


£ 


d 


-55.0 


X 


0.8 


-150.3 


± 


1.9 


-8.37 


X 


0.22 


53.0 


X G 


T A c* r* 
1 Aul 


1 A C 




c 


-58.7 


X 


2.1 


-162.3 


x 


6.5 


-8.39 


X 


0.10 


£1 A 

52.0 


cc 


ATCG 


A T G 


1 


d 


-57.4 


± 


1.9 


-158.2 


± 


4.6 


-8.32 


X 


0.48 


52.0 


£ G 


T A G C 


T A C 


I 




-59.3 


db 


1.6 


-160.7 


± 


5.1 


-8.21 


X 


0.07 


58.1 


1 C 


ATCG 


AT G 


£ 


d 


-59.5 


db 


1.2 


-165.4 


x 


2.3 


-8.19 


X 


0.28 


50.7 




n A T A 

yj A I A 


i-p /-i 


£ 


c 


-62.8 


± 


1.2 


-175.7 


x 


3.9 


-8.34 


X 


0.12 


50.7 


£ G 


C T A T 


AG C 


I 


d 


7 

O /. / 


x 






L 
X 


J. / 


-o. 10 


) 

a; 




51.0 




A T A 


t r* c 


X 


c 


-61.1 


± 


1.8 


-170.6 


X 


5.5 


-8.13 


± 


0.13 


50.0 


I G 


C T A T 


A G C 


Q 


d 


-56.3 


± 


1.3 


-155.8 


X 


3.1 


-7.98 


± 


0.34 


50.2 


I A 


GAGC 


T C T 


c 


c 


-56.4 


± 


2.0 


-153.5 


X 


6.3 


-8.81 


± 


0.12 




£ T 


CTCG 


A G A 




d 


-53.8 


± 


1.8 


-145.5 


X 


4.3 


-8.69 


± 


0.52 


55.4 


£ A 


GAGC 


T C T 


i 


c 


-54.8 


± 


2.1 


-149.4 


X 


6.4 


-8.47 


db 


0.17 


53.7 


I T 


CTCG 


A G A 


£ 


d 


-50.3 


db 


1.2 


-135.6 


X 


2.8 


-8.28 


db 


0.36 


53.9 
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Tabl&Sl: Continued. 





























AH° 


AS' 


0 


AG° 37 








(kcal / mol) 


(cal/molK) 


(kcal / mol) 


(°C) 


fiTGAGCTCAA 


c 


-56.6 ± 


2.8 


-154.3 


± 


8.4 


-8.79 




0.17 


55.1 


AACGCG AGT£ 


d 


-53.5 ± 


1.6 


-144.6 


± 


3.6 


-8.64 


± 


0.44 


55.2 


ATGAGCTCAG 


c 


-57.7 ± 


1.9 


-157.0 


± 


5.6 


-9.05 


± 


0.15 


56.3 


SACTCG AGT A 


d 


-53.8 ± 


1.6 


-145.1 


± 


3.6 


-8.85 


± 


0.45 


56.5 


fiGTAGC TACA 


c 


-59.2 ± 


4.2 


-162.7 


a. 


12.9 


-8.77 




0.16 


54.1 


AC ATCGATGG 


d 


-59.4 ± 


1.7 


-163.3 


± 


4.1 


-8.72 




0.44 


53.8 


AGTAGCTACG 


c 


-63.3 ± 


3.4 


-173.4 


± 


10.2 


-9.46 


± 


0.20 


56.8 


GCATCGATGA 


a 


-63.8 ± 


1.8 


-175.2 


± 


4.5 


-9.47 


± 


0.45 


56.6 


fiCGATATCGA 


c 


-62.1 ± 


1.3 


-172.8 




3.8 


-8.51 


± 


0.10 


51.8 


AGCTAT AGC£ 


a 


-59.1 ± 


1.1 


-163.6 




.2.6 


-8.39 




0.27 


51.9 


ACGATATCGfi 


c 


-63.1 ± 


2.6 


-174.1 


± 


7.9 


-9.06 


± 


0.15 


54.6 


GGCTAT AGC A 


d 


-60.4 ± 


1.6 


-166.0 




3.7 


-8.92 




0.40 


54.6 


SAGAGCTCTA 


c 


-56.4 ± 


1.9 


-153.5 


± 


5.6 


-8.84 




0.15 


55.4 


AT CT CG AGA£ 


d 


-52.8 ± 


1.3 


-142.4 


± 


3.1 


-8.67 


± 


0.38 


55.6 


AAGAGCTCTfi 


c 


-57.1 ± 


2.6 


-154.7 


± 


7.6 


-9.08 




0.20 


56.7 


GT CTCGAGA A 


d 


-53.1 ± 


2.2 


-142.8 


db 


6.7 


-8.86 


± 


0.10 


56.8 


IT GAGC T CA£ 


c 


-55.2 ± 


1.5 


-150.0 


± 


4.5 


-8.63 


± 


0.07 


54.5 


fiACTCG AGT 1 


d 


-57.2 ± 


0.8 


-156.4 




1.8 


-8.71 


± 


0.21 


54.3 


GTGAGCTCAI 


c 


-57.3 ± 


0.7 


-156.5 


± 


2.3 


-8.80 


± 


0.05 


54.9 


lACTCGAGTfi 


A 
U 


-57.6 ± 


1.6 


-157.4 


± 


3.7 


-8.81 


± 


0.42 


54.8 


IGTAGCTACfi 


c 


-59.9 ± 


1.8 


-164.5 


± 


5.5 


-8.93 


± 


0.08 


54.8 




d 


-58.2 ± 


2.0 


-158.9 


± 


4.8 


-8.86 


± 


0.53 


55.0 


fiGTAGC T AC 1 


c 


-62.8 ± 


3.6 


-174.8 


± 


11.0 


-8.62 




0.17 


52.3 


ICATCGATGfi 


d 


-67.1 ± 


3.8 


-187.8 


± 


9.4 


-8.79 


± 


0.85 


52.1 
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Table SI: Continued. 





AH 0 
(kcal / mol) 


(cal/molK) 


AG 0 ™ 
(kcal / mol) 


(°C) 


TTGATATCGG 
fiCCTAT AGC1 


c -58.9 ± 
d -63.5 ± 


2.8 
0.6 


-163.0 ± 8.9 
-177.3 ± 1.6 


-8.40 ± 0.06 
-8.53 ± 0.15 


52.0 
51.6 


GCGATATCGT 
1GCTAT AGCfi 


c -61.6 ± 


3.7 


-170.4 ±11.5 
-177.6 ± 2.3 


-8.71 ± 0.14 
-8.77 ± 0.22 


53.1 
52.8 


XAGAGCTCTfi 
fiTCTCGAGAl 


c -56.4 ± 
d -58.0 ± 


1.5 
0.6 


-153.3 ± 4.5 
-158.1 ± 1.9 


-8.89 ± 0.06 
-8.96 ± 0.03 


55.7 

- - e 


Q A GAGCTCTI 
XT CTCGAGAS 


c -57.7 ± 
d -56.6 ± 


1.6 
1.9 


-157.2 ± 4.7 
-153.9 ± 4.6 


-8.95 ±0.10 
-8.90 ± 0.52 


55.7 
55.7 


Core sequences 

C GATATCG' 
GCTAT AGC 


c -55.7 ± 
d -51.8 ± 


3.9 
0.6 


-157.1 ± 12.1 
-145.1 ± 1.4 


-6.93 ± 0.12 
-6.82 ± 0.15 


44.1 
44.0 


GTAGC T AC' 
CATCGATG 


c -55.1 ± 
d -51.4 ± 


2.3 
0.6 


-155.0 ± 7.0 
-143.3 ± 1.3 


-7.04 ± 0.10 
-6.95 ± 0.14 


44.9 
44.9 


AGAGCTCT 
T CTCGAGA 


c -49.5 -± 
d -50.1 ± 


1.8 
0.8 


-134.5 ± 5.7 
-136.7 ± 1.8 


-7.76 ± 0.07 
-7.76 ± 0.22 


50.6 
50.5 


T GAGCTCA 
ACTCGAGT 


c -50.7 ± 
d -50.3 ± 


0.7 
0.7 


-138.4 ± 2.2 
-137.3 ± 1.6 


-7.73 ± 0.04 
-7.72 ± 0.18 


50.1 
50.1 



a The top strand of each duplex is represented in the 5' to 3' orientation and the bottom strand 
is shown in the 3' to 5' direction. Terminal mismatch jiearest neighbors are represented in bold. 
Mismatches are underlined. b T M calculated using 10" 4 M total strand concentration. 

c Thermodynamic parameters from averaging the fits of melting curves. Reported errors are 
standard deviations in the precision of the data. a Thermodynamic parameters from T M * vs. in(C T ) 

plots. Reported errors are standard deviations in the precision propagated from the slope and 
intercept of the 1/T M vs. In C T plot. c Data from reference (19). 
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Table 1 : Thermodynamic Parameters for Hairpin Oligomer Association and Oligomer Dupiex 
Formation. 

AH° AS 0 AG° r t m 

(kcal/mol) (caI/'mo!K) (kcaJ mol) i°Ci 

Systems with Elementary Interfaces'* 

-AAGCCTTGA-ACAACG c -50.7 ± 4.1 -141.7 ± 11.3 -6.74 = 0.27 49.3 
KGCGGAACT/TGTTGC (i) 



r AAGCCTTGT-TCAACG c -52.8 ± 4.2 -146.3 * 11.7 -7.43 = 0.30 53.1 
HGCGGAACA/AGTTGC 



GCAACT-TGTTCCGAAn e -63.5 ± 5.1 -179.2 ± 14.3 -7.99 ± 0.32 53.2 
C G TTG A / A C AAGGCCCV 

-AAGCCTT GA - TCAACG e -53.6 ± 4.3 -149.3 ± 11.9 -7.34 i 0.29 52.3 
^-CGCGGAACT/AGTTGC 

AAGCCTT GT - ACAACG e -45.1 ± 3.6 -124.8 ± 10.0 -6.42 ± 0.26 48.4 
CGCGGAACA/TGTTGC 

AAGCCTT GC -ACAACG e -46.1 ± 3.7 -128.5 ± 10.3 -6.26 * 0.25 46.9 
CGCGGAACG/TGTTGC 

AAGCCTT GT - GCAACG c -52.2 ± 4.2 -144.4 ± 11.5 -7.39 ± 0.30 53.1 
CGCGGAACA/CGTTGC 

A AGC C TT GG - TCAACG c -53.6 ± 4.3 -148.2 ± 11.9 -7.67 * 0.31 54.4 



C 



C 



C 



C 



CGCGGAACC / AGTTGC 



/-AAG C C T T G A - CC AA C G e -51.3 * 4.1 -140.2 * 11.2 -7.81 i 0.31 56.2 
^■CGCGGAAC T /GGTTGC 
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Table 1 : Continued. 



AH 0 AS 0 AGV t m 

(kcal/mol) (cal/molK) (kcal moh ( °C> 

c A AGC C T T G C - TC AAC G e -46.1 ± 3.7 -126.3 ± 10.1 -6.90 = 0.28 51." 
^CGCGGAA CG / AGTTGC 



,-AAGCCTT GA-GCAACG e -48.2 ± 3.9 -131.7 ± 10.5 -7.32 ± 0.29 53.9 
^CGCGGAACT/CGTTGC 



GCAACA «G GTTCC GAA-v e -51.9 ± 42 -147.0 ± 11.8 -6.34 ± 0.25 46.3 
CGTTGT / C CAAGGCC C-' 



.AAGCCTT GG - ACAACG e -50.4 ± 4.0 -139.5 ± 11.2 -7.10 ± 0.28 51.7 
GCGGAA C C / TGTTG C 



/•AAGCCTT GT - CC AAC G c -54.2 ± 4.3 -147.9 ± 11.8 -8.29 ± 0.33 58.3 
^CGCGGAA CA/GGTTGC 



/-AAGCCTT GC-GCAACG c -47.6 ± 3.8 -130.9 ± 10.5 -6.96 ± 0.28 51.6 
^CGCGGAA CG/CGTTGC 



.-AAGCCTT GG-CCAACG e -52.1 ±4.2 -140.1 ± 11.2 -8.67 ± 0.35 61.9 
^•CGCGGAA CC/GGTTGC 



r A AGC CTT GC - CCAACG e -53.3 ± 4.3 -146.5 ± 11.7 -7.92 ± 0.32 56.1 
^•CGCGGAA CG/GGTTGC 



--AAGCCTT GG-GCAACG e -49.3 ± 3.9 -1352 ± 10.8 -7.35 ± 0.29 53.8 
^■CGCG G A ACC/CGTTGC 
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Table 1: Continued. 



AH° 
(kcal / mol) 

Systems with Dangling Ends at th>» Tnt»rfan. b 

.AAGCCTTGC-GCAACG c -44.4 ± 3.6 
^CGCGGAA CG / CGTTGC 
A 

.AAGCCTTGC-GCAACG c -48.0 ± 3.8 
^CGCGGAACG/CGTTGC 

A 

.AAGCCTTGCGCAACG c -46.3 ± 3.7 
GCGGAA C G / CGTTGC 
T 

^AAGCCTT GC-GCAACG c -49.0 ± 3.9 
^•C GCGGAA C G / CGTTGC 

T 

,AAGCCTT GC-GCAACG e -37.6 ± 3.0 
^-C GCGGAA C G / CGTTGC 
A A 



(cal / mol K) (kcal / mol) i*c ) 



-121.4 ± 9.7 -6.76 ± 0.27 51.1 



-131.8 ± 10.5 -7.16 ± 0.29 52.9 



-127.4 ± 10.2 -6.83 ± 0.27 51.0 



-136.6 ± 10.9 -6.59 ± 0.26 48.7 



•102.0 ± 8.2 -5.91 ± 0.24 46.4 



-AAGCCTTGC - GCAAC G c -36.2 ± 2.9 -97.3 ± 7.8 -6.03 ± 0.24 47.8 
•C GCGGAA C G / CGTTGC 
T T 

-A AGCCTT GC-GCAACG e -44.0 ± 3.5 -123.5 ± 9.9 -5.67± 0.23 43.0 
•C GCGGAA CG / CGTTGC 
A T 

• A AGCCTT GC- GCAAC G c -43.6 ± 3:5 -119.5 ± 9.6 -6.53 ± 0.26 49.6 
•CGCGGAACG/CGTTGC 
T A 

.AAGCCTTGG-TCAACG e -47.2 ± 3.8 -129.8 * 10.4 -6.90 ± 0.28 51.3 
•CGCGGAACC/AGTTGC 

A 
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Table 1: Continued. 



AH 0 AS 0 AG° 37 T M 

(kcal/mol) (cal/molK) (kcal/mol) ( °C) 

Svstems with Extra Central Nucleotide at t he Interface" 

^AAGCCTT GCAGCAACG c -44.4 ± 3.6 -1222 ± 9.8 -6.50 * 0.26 49.2 
GCGGAA C G / CGTTGC 

r A AGC CTT GTACCAACG 6 -45.0 ± 3.6 -124.6 ± 10.0 -6.32 ± 0.25 47.7 
^CGCGGAACA/GGTTGC 



Oligomers 

TCAACG c -38.5 ± 3.1 -108.3 ± 8.7 -4.94 ± 0.20 38.0 
AGTTGC 



AC AA C G c -36.1 ± 2.9 -99.3 ± 7.9 -5J2 ± 0.21 41.4 
TGTTGC 



GCAACG c -42.5 ± 3.4 -117.1 ± 9.4 -621 ± 025 47.5 
CGTTGC 



CCAACG c -38.6 ± 3.1 -106.2 ± 8.5 -5.69 ± 023 44.2 
GGTTGC 



TGTTGC c -37.1 ± 3.0 -101.1 ± 8.1 -5.79 ± 023 45.3 
ACAACG 



AGTTGC e -37.0 ±3.0 -101.0 ± 8.1 -5.70 ± 023 44.2 
TCAACG 
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a T M calculated using 10" 4 total strand concentration. 

b The top strand of each system is conventionally represented in the 5' to 3' orientation. Nucleotides 
involved in coaxial stacking interfaces are represented in bold. 

c Parameters obtained by averaging the results of melt fit and TM" 1 vs. ln(C T /4) plot methods. 
Errors are estimated to be 8% for AH 0 and AS 0 and 4% for AG° 37 . 
(0. (id. (no Melt j ng curV es for these systems are shown in Figure 2. 
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Table 2: Thermodynamic Parameters for Coaxial Stacking*. 





AH°(coaxial stacking) 
(kcal / mol) 


AS°(coaxial stacking) 
(cal / mol K) 


AG° 37 ( coaxial stacking) 
(kcal / mol) 


Klementarv Interfaced 






GA - AC 
CT / TG 


-14.6 ± 5.0 


-42.4 ± 13.8 


-1.42 ± 0.34 


GT-TC 
CA / AG 


-14.3 ± 5.2 


-38.0 ± 14.6 


-2.49 ± 0.36 


CT - TG 
GA / AC 


-26.6 ± 5.9 


-78.2 ± 16.5 


-2^9 ± 0.39 


GA - TC 
CT / AG 


-15.1 ± 53 


-41.0 ± 14.8 


-2.40 ± 035 


GT - AC 
CA / TG 


-9.0 ± 4.6 


-25.5 ± 12.8 


-1.10 ± 033 


GC - AC 
CG / TG 


-10.0 ± 4.7 


-29.2 ± 13.0 


-0.94 ± 033 


GT - GC 
CA / CG 


-9.6 ± 5.4 


-27.3 ± 14.9 


-1.18 ± 039 


GG - TC 
CC / AG 


-15.1 ± 5.3 


-39.9 ± 14.7 


-2.73 ± 0.36 


GA - CC 
CT / GG 


-12.7 ± 5.1 


-34.0 ± 14.1 


-2.12 ± 0.39 


GC - TC 
CG / AG 


-7.6 ± 4.8 


-18.0 ± 133 


-1.97 ± 034 


GA-GC 
CT / CO 


-5.6 ± 5.1 


•14.6 ± 14.1 


-1.11 ± 0.38 


CA - GG 
GT / CC 


-14.8 ± 5.1 


-45.9 ± 143 


-0.56 ± 034 


GG - AC 
CC / TG 


-14.2 ± 5.0 


^0.2 ± 13.7 


-1.78 ± 0.35 
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Table 2: Continued. 



GT - CC 
CA/GG 

GC - GC 
CG / CG 

GG - CC 
CC / GG 

GC - CC 
CG / GG 

GG - GC 
CC / CG 



AH°(coaxiai stacking) 
(kcal/mol) 

-15.6 ± 5.3 



-5.0 ± 5.1 



•13.5 ± 52 



-14.7 ± 5.3 



-6.8 ± 52 



AS°(coaxial stacking) AG°3 7 (coaxial stacking) 
(cal/molK) (kcal/mol) 



Interfaces with Panning Ends b 



GC - GC 
CG / CG 
A 

GC - GC 
CG / CG 
A 

GC - GC 
CG / CG 
T 

GC - GC 
CG / CG 
T 

GC - GC 
CG / CG 
A A 



■1.9 ± 4.9 



-5.5 ± 5.1 



-3.8 ± 5.0 



-6.4 ± 52 



5.0 ± 4.5 



-41.8 ± 14.6 



•13.8 ± 14.1 



-33.9 ± 14.1 



-40J ± 14.5 



•18.1 ± 14.3 



-4.3 ± 13.5 



-14.7 ± 14.1 



-10J ± 13.8 



-19.5 ± 14.4 



15.1 ± 12.4 



-2.61 ± 0.40 



-0.75 ± 0.37 



-2.98 ± 0.41 



-2.23 ± 0.39 



-1.14 i 0.38 



-0.55 ± 0.37 



-0.95 ± 0.38 



-0.62 ± 0.37 



-0.38 ± 0.36 



0.30 ± 0.34 
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Table 2: Continued. 



AH°(coaxial stacking) AS°(coaxial stacking) AG° 37 ( coaxial stacking) 

(kcal/mol> (cal/molK) (kcal/tnol) 

GC-GC 6.3 ± 4.5 19.8 ± 12.2 0.18 ± 0.35 

CG / CG 
T T 

GC-GC -1.4 .± 4.9 -6.4 ± 13.6 0.55 ± 0.34 

CG / CG 
A T 

GC-GC -1.1 ± 4.9 -2.4 ± 13.4 -0.32 ± 0.36 

CG/CG 

T A 

A 

GG-TC -8.6 ± 4.9 -21.5 ± 13.5 -1.96 ± 0.34 

CC/AG 
A 

Interface with Extra Central Nuelentirie* 

GCAGC -1.9 ± 4.9 -5.1 ± 13.5 -0.29 ± 0.36 

CG/CG 

GTACC -6.4 ± 4.7 -18.5 ± 13.1 -0.64 ± 0.34 

CA / GG 



* These parameters and their corresponding errors are deduced from Table 1 as described in the text 
b The top strand of each duplex is conventionally represented in the 5' to 3' orientation. Nucleotides 
involved in coaxial stacking interfaces are represented in bold. 
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Table SI: Extinction coefficients of hairpins at 25 °C 



C 



AAGCCTTGGTCAACG 
C GCGGAACC 



evpgrimental* calculated* 
(L mol" 1 cm" 1 ) (L mol" 1 cm*') 



188847 



192310 



C 



AAGCCTTGATCAACG 
C GCGGA ACT 



188718 



195950 



C 



AAGCCTTGGTCAACG 
CGCGGAACG 



186139 



191610 



C 



AAGCCTTGTTCAACG 
C GCGGA ACA 



191071 



195810 



AAGCCT TGAACAACG 
■C GCGGA ACT 



194330 



200950 



r AAGCCTTGTACAACG 
^CGCGGAACA 



193953 



200810 



r A AGCC T TGGACAACG 
^-C GCGGAACC 



188889 



197310 



AAGCCTTGCACAACG 
CGCGGAACG 



194623 



196610 



c 



A AGCC T TGAG.CAACG 
C GCGGAACT 



192953 



197350 



c 



AAGCCTTGTGC AAC G 
C GCGGAACA 



195968 



197630 



A AGCC T TGGGCAACG 
C GCGGAACC 



195177 



193710 
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Table SI: Continued 



C 



A A G C C T T GCGCAAC G 
CGCGGAACG 



190944 



193010 



c 



AAGCCTTGACCAACG 
C GCGGAACT 



192663 



194350 



c 



A AGC C TTGTCCAACG 
CGCGGAACA 



193532 



194210 



AAGC CTTGGCCAACG 
•C GCGGAACC 



195094 



192390 



C 



AAGC CTTGCCCAACG 
CGCGGAACG 



192864 



190010 



GCAACA-G TTCC A A 
CCAAGGCCC 



200806 



196010 



GCAACT-T TTCC AAn 
ACAAGGCcJ 



206650 



192810 



•AAGCCTTGCAGCAACG 
•CGCGGAACG 



202688 



206510 



C 



AAGCCTTGTACCAACG 
CGCGGAACA 



204450 



208010 



AAGC CTTGCGCAACG 
CGCGGAACG 
A 



191282 



193010 



AAGCCTTGCGCAACG 
CGCGGAACG 
T 



198343 



193010 
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3 Calculated with Equation 3? 
b Calculated with Equation 4. 
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Table S2: Thermodynamic Parameters for Hairpin Oligomer Association and Oligomer Duplex 
Formation. 



AH 0 AS 0 AG° 3: T M 3 

(kcal/mol) (cal/mol.K) (kcal/mo!) ( °C) 

pigmentary interfaces' 



AAGCCTTGA-ACAACGW « 
GCGGAACT / TGTTGC c 


-48.1 ± 
1 -53.3 ± 


2.5 
2.7 


-133.2 i 
-150.1 ± 


8.1 
8.7 


-6.8 
-6.7 


± 


0.1 
0.0 


50.1 
48.4 


-AAGCCTTGT - T C A A C G <»> « 
^CGCGGAACA/AGTTGC « 


■ -55.7 ± 
1 -49.9 ± 


4.3 
1.7 


-155.8 ± 
-136.9 ± 


14.0 
5.6 


-7.4 
-7.4 


± 
± 


0.2 
0.0 


52.1 
54.2 


r AAGCCTTGA-TCAACG (TO < 
^CGCGGAACT / AGTTGC 


» -51.9 ± 
1 -55.4 ± 


4.0 
2.1 


-143.6 ± 
-155.0 ± 


12.7 
4.9 


-7.3 
-7.3 


± 
± 


0.1 
0.5 


52.8 
51.7 


GCAACT-TGTTCCGAA-v « 
CGTTGA / ACAAGGCCC-/ * 


' -61.6 ± 
1 -65:5 ± 


2.6 
3.1 


-173.1 ± 
-185.2 ± 


7.9 
9.8 


-8.0 
•8.0 


± 
± 


0.1 
0.0 


53.7 
52.8 


-AAGC C TTGT - AC A A C G « 
^CGCGGAACA/ TGTTGC ' 


■ -48.0 ± 
1 -42.2 ± 


2.3 
1.9 


-134.4 ± 
-115.2 ± 


7.2 
6.3 


-6.4 
-6.5 


± 
± 


0.1 
0.1 


47.3 
49.5 


-AAGCCTTGC-ACAACG « 
<-C GCGGAACG / TGTTGC * 


» -45.9 ± 
1 -46.3 ± 


1.1 
4.0 


-127.9 ± 
-129.2 ± 


3.3 
13.1 


-6.3 
-6.3 


± 
± 


0.1 
0.1 


47.1 
46.8 


r AAGCCTTGT-GCAACG ' 
^CGCGGAACA/ CGTTGC ' 


• -54.7 ± 
1 -49.7 ± 


2.7 
3.4 


-152.4 ± 
-136.3 ± 


9.3 
11.0 


-7.4 
-7.4 


± 
± 


0.2 
0.1 


52.2 
53.9 


r AAGCCTTGG- TCAACG « 
^C GCGGAACC / AGTTGC ' 


» -53.5 ± 
1 -53.8 ± 


4J 
3.4 


-147.7 ± 
-148.7 ± 


13.6 
11.1 


-7.7 
-7.7 


± 
± 


0.2 
0.1 


54.4 
54.3 


r AAGCCTTGA- CCAACG « 
^-CGCGGAACT / GGTTGC ' 


' -52.9 ± 
1 -49.6 ± 


2.2 
2.4 


•145.4 ± 
-134.9 ± 


7.3 
5.5 


-7.8 
-7.8 


± 
± 


0.1 
0.7 


55.6 
56.8 


r AAGCCTTGC-TCAACG < 
^C GCGGAACG / AGTTGC ' 


e -48.6 ± 
1 -43.6 ± 


2.8 

2.2 


-134.6 ± 
-118.1 ± 


9.3 
7.3 


-6.9 
-6.9 


± 
± 


0.2 
0.1 


50.6 
52.8 


r AAGCCTTGA-GCAACG 
^CGCGGAACT/ CGTTGC 


' -49.7 ± 
4 -46.6 ± 


3.5 
1.8 


-136.6 ± 
-128.8 ± 


11.4 
5.7 


-7.3 
-7.3 


± 
± 


0.1 
0.0 


53.4 
49.9 
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Table S2: Continued. 



GCAACA - GGTTCCGA A-v 
CGTTGT / CCAAGGCCC-' 

-AAGCCTTGG-ACAACG 
^-C GCGGAACC / TGTTGC 

-AAGCCTTGT - CCAACG 
^CGCGGAACA/GGTTGC 

-AAGCCTTGC-GCAACG 
^•CGCGGAACG / CGTTGC 

-AAGCCTTGG-CCAACG 
^CGCGGAACC/GGTTGC 

r AAGCCTTGC- CCAACG 
^»C GCGGAACC / GGTTGC 

-AAGCCTTGG-GCAACG 
^C GCGGAACC /CGTTGC 



AH 0 
(kcal / mol) 

-53.1 ± 3.8 
-50.8 ± 3.3 

-48.4 ± 3.1 
-52.3 ± 2.0 

-58.0 ± 2.1 
-50.3 i 2.1 

-47.0 ± 3.9 
-48.2 ± 0.9 

-57.3 ± 5.4 
-46.9 ± 2.7 



-55.3 ± 
-51.4 ± 

-45.9 ± 
-52.6 ± 



2.6 
1.5 

1.0 
3.5 



AS 0 
(cal / mol K) 

-150.8 = 13.0 
•143.3 ± 10.9 



-133.1 
-145.9 

-160.2 
-135.7 



-129.0 ± 

-132.8 ± 

•156.4 ± 

-123.9 ± 

•152.7 ± 

-140.2 ± 

-124.4 ± 

-146.0 ± 



10.4 
6.5 

6.4 
4.7 

12.9 
2.9 

17.0 
5.8 

8.1 
3.4 

3.5 
113 



AGV V 
(kcal -mol) t °Ci 



-6.3 = 0.2 

-6.4 = 0.1 

-7.1 £ 0.1 

-7.1 i 0.0 

-8.4 ± 0.2 

-8.2 ± 0.6 

-7.0 ± 0.1 

-7.0 ± 0.0 



45.8 
46.9 

52.5 
50.9 

57.1 
59.5 

51.7 
51.5 



-8.8 ± 0.3 60.1 
-8.5 ± 0.9 63.6 



-7.9 
-7.9 



0.1 
0.5 



-73 ± 0.1 
-7.4 ± 0.1 



55.6 
56.7 

54.9 
52.7 



Interfaces with dangling ends* 



.AAGC CTTGC - GCAACG 
CGCGGAACG/ CGTTGC 
A 



c 



AAGC CTTGC - GCAACG 
CGCGGAACG/ CGTTGC 
T 

AAGCCTTGC-GCAACG 
CGCGGAACG/ CGTTGC 

A 



-453 ± 2.9 
-43.7 ± 3.4 



-46.7 ± 3.7 
-46.0 ± 33 



-48.6 ± 3.0 
-47.4 ± 33 



-123.8 ± 9.4 
-119.1 ± 11.1 



•128.7 ± 11.8 
-1263 ± 7.5 



•133.7 ± 9.4 
•129.8 ± 10.5 



-6.8 ± 03 50.8 
-6.8 ± 0.1 51.5 



-6.8 ± 0.1 50.9 
-6.8 ± 1.0 513 



-73 ± 0.1 52.7 
-73 ± 0.1 53.1 
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Table S2: Continued. 



-A AGCCTTGC- GCAACG 
^-C GCGGAACG / CGTTGC 



c 



AAGCCTTGG- TCAACG 
CGCGGAACC/AGTTGC 

A 



r A AGC C TTGC-GCAACG 
S C GCGGAACG/CGTTGC 
A A 



r A AGC CTT GC - GCAAC G 
^CGCGGAACG/CGTTGC 
A T 



A AGCCTTGC- GCAACG 
CGCGGAACG/CGTTGC 
T T 



r AAGCCTTGC-GCAACG 
^CGCGGAACG/CGTTGC 
T A 

Interfaces with extra central nucleotide 



AH° 
(kcal / mol) 

-46.3 ± 2.6 
-51.6 ± 3.7 



AS° 
(cal / mol K) 

-127.9 i 8.7 
-145.4 ± 12.1 



AGV V 
(kcal /mol) j°C) 



-6.6 * 0.2 
-6.5 = 0.1 



43.6 
42.3 



-49.2 ± 3.5 -136.6 ± 11.0 -6.9 ± 0.2 44.8 
-45.1 ± 2.9 -123.1 ± 9.4 -6.9 ± 0.1 45.7 



-39.7 ± 
-35.5 ± 



-44.7 ± 
-432 ± 



-38.5 ± 
-33.9 ± 



4.6 
2.9 



33 
3.7 



3.0 
13 



-43.3 ± 9.0 
-43.9 ± 1.0 



-109.1 ± 
-95.0 ± 



-126.0 ± 
-121.0 ± 



-105.1 
-89.6 



-118.5 
-120.5 



15.2 
5.9 



11.0 
12.2 



10.0 
4.5 



30.1 
3.4 



-5.8 ± 0.2 38.2 
-6.0 i 1.1 39.8 



-5.6 
-5.7 



-6.0 
-6.1 



-6.5 ± 
-6.5 ± 



0.2 
0.2 



0.2 
0.1 



0.3 
0.0 



36.6 
37.2 



39.3 
41.0 



43.1 
432 



r A AGC CTTGCAGCAACG 
^•CGCGGAACG/CGTTGC 



« -43.8 ± 9.4 
4 -45.0 ± 3.5 



-1202 ± 
-124.1 ± 



31.0 
11.5 



-6.5 ± 
-6.5 i 



0.2 
0.1 



43.0 
42.7 



AAGCCTTGTACCAACG « -45.0 ± 2.8 



CGCGGAACA/GGTTGC 



-45.0 ± 4.0 



-124.5 
-124.8 



9.1 
13.4 



-63 ± 0.2 
-63 ± 0.2 



41.6 
41.4 
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Table S2: Continued. 



AH 0 AS 0 AG 0 ,, T M J 

(kcal/mol) (cal/molK) (kcal.'mol) { *C\ 

QligQmfiB 



T C A A C G 


e -40.6 ± 


2.5 


-115.4 


± 


9.1 


-4.8 




0.3 


30.5 


AGTTGC 


d -36.5 ± 


1.9 


-101.2 


± 


6.3 


-5.1 


± 


0.1 


31.9 


. At A A L Cr 


c -3o.J X 


z.u 


1 ML o 

•lUo.o 


* 

± 


< < 

0.3 






U.- 




TGTTdC 


d -33.9 ± 


1.6 


-91.8 


S 


5.3 


-5.4 


X 


O.I 


34.6 • 


GC AAC G 


c -43.9 ± 


2.4 


-121.6 


X 


7.4 


-6.2 


± 


0.1 


40.7 


CGTTGC 


d -41.1 ± 


2.1 


-112.6 


X 


6.9 


-6.2 


X 


0.1 


41.2 


CCAACG 


e -41.0 ± 


2.6 


-114.1 


X 


8.9 


-5.6 


X 


0.2 


36.4 


GGTTGC 


d -362 ± 


0.9 


-98.2 


X 


1.8 


-5.8 


X 


0.3 


37.8 


TGTTGC 


c -38.8 ± 


3.8 


-106.6 


X 


12.1 


-5.7 


X 


0.1 


37.5 


A C A A C G 


d -35.5 ± 


2.8 


-95.7 


X 


9.2 


-5.8 


X 


0.1 


38.4 


AGTTGC 


c -37.1 ± 


2.9 


-101.2 


± 


9.4 


-5.7 


X 


0.1 


36.9 


TCAACG 


d -36.9 ± 


2.4 


-100.8 


± 


8.0 


-5.7 


X 


0.1 


36.9 



a T M calculated for 4x1 0* 4 total strand concentration. 

b The top strand of each system is conventionally represented in the 5* to 3' orientation. Nucleotides 
involved in coaxial stacking interfaces are represented in bold. 

c Parameters obtained from averaging fits of melting curves. Reported errors are standard deviations 

in the precision of the fitted data. 
Parameters obtained from T M vs. ln(CV4) plots. Reported errors are standard deviations in the 

precision propagated from the slope and intercept of the 1/T M vs. In (<V4) plot 

OKffKOQ iyr M vs. to(C|/4) plots for these systems are shown in Figure SI. 
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WHAT IS CLAIMED IS: 



1 LA method for predicting nucleic acid hybridization 

2 thermodynamics, the method comprising: 

3 providing a database of thermodynamics parameters; 

4 receiving hybridization information which represents at least one 

5 sequence; 

6 receiving correction data; 

7 receiving a first set of data which represents hybridization conditions; 

8 and 

9 calculating hybridization thermodynamics including net hybridization 

10 thermodynamics based on the hybridization information, the thermodynamic 

1 1 parameters, the correction data and the first set of data. 

1 2. The method as claimed in claim 1 wherein the hybridization 

2 thermodynamics of individual single stranded, bimolecular and higher order 

3 complexes are statistically weighted in a numerical process and the equilibrium 

4 concentration of each species is output. 

1 3 . The method as claimed in claim 2 wherein the correction data 

2 includes folding correction data. 

1 4. The method as claimed in claim 2 wherein the correction data 

2 includes linear correction data. 

1 5 . The method as claimed in claim 1 wherein the thermodynamic 

2 parameters include DNA thermodynamic parameters. 

1 6. The method as claimed in claim 5 wherein the DNA 

2 thermodynamic parameters include dangling end parameters. 

1 7. The method as claimed in claim 5 wherein the DNA 

2 thermodynamic parameters include coaxial stacking parameters. 
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1 8. The method as claimed in claim 5 wherein the DNA 

2 thermodynamic parameters include terminal mismatch parameters. 

1 9 . The method as claimed in claim 1 wherein the thermodynamic 

2 parameters include RNA thermodynamic parameters. 

1 10. The method as claimed in claim 1 wherein the thermodynamic 

2 parameters include hybrid DNA/RNA thermodynamic parameters. 

1 11. The method as claimed in claim 1 wherein the thermodynamic 

2 parameters include DNA loop thermodynamic parameters. 

1 12. The method as claimed in claim 1 wherein the hybridization 

2 information represents top and bottom strand sequences which form a duplex and 

3 wherein the hybridization thermodynamics are calculated for the duplex. 

1 13. The method as claimed in claim 1 wherein the hybridization 

2 information represents at least a section of a target and a length of at least one 

3 primer or probe complimentary to the target. 

1 14. The method as claimed in claim 13 wherein the hybridization 

2 thermodynamics are calculated for a plurality of primers or probes complimentary 

3 to the target. 

1 15. The method as claimed in claim 1 wherein the hybridization 

2 information represents at least a section of a target and a primer or probe. 

1 16. The method as claimed in claim 15 wherein a length of the 

2 target is longer than a length of the primer or probe and wherein the hybridization 

3 thermodynamics are calculated for a best target/primer or target/probe complex and 

4 for competitive mismatch complexes. 
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1 17. The method as claimed in claim 14 wherein hybridization 

2 information represents at least a section of a target and a primer or probe and 

3 wherein a length of a target is longer than the length of the primer or probe and 

4 wherein the hybridization thermodynamics are calculated for a best target/primer 

5 or target/probe complex and for competitive target/primer or target/probe 

6 complexes. 

1 18. The method as claimed in claim 2 further comprising, 

2 calculating concentration of each species in a solution at a plurality of temperatures. 

1 19. The method as claimed in claim 18 wherein hybridization 

2 information also represents a primer or probe and wherein the length of the target 

3 is longer than a length of the primer or probe and wherein the hybridization 

4 thermodynamics are calculated for a best target/primer or target/probe complex arid 

5 ; for competitive mismatch complexes and wherein the method further comprises 

6 calculating concentration of every species in a solution at a plurality of 

7 temperatures. 

1 20. The method as claimed in claim 19 wherein the hybridization 

2 thermodynamics are calculated for at least two best target/primer or target/probe 

3 complexes and for their corresponding competitive mismatch complexes and 

4 wherein the method further comprises correcting for any interactions between the 

5 at least two best target/primer or target/probe complexes and their components. 

1 21. A system for predicting nucleic acid hybridization 

2 thermodynamics, the system comprising: 

3 a database of thermodynamics parameters; 

4 means for receiving hybridization information which represents at 

5 least one sequence; 

6 means for receiving correction data; 

7 receiving a first set of data which represents hybridization conditions; 

8 and 
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9 means for calculating hybridization thermodynamics including net 

10 hybridization thermodynamics based on the hybridization information, the 

1 1 thermodynamic parameters, the correction data and the first set of data. 

1 22. The system as claimed in claim 21 wherein the hybridization 

2 thermodynamics of individual single stranded, bimolecular and higher order 

3 complexes are statistically weighted in a numerical process and the equilibrium 

4 concentration of each species is output. 

1 23 . The system as claimed in claim 22 wherein the correction data 

2 includes folding correction data. 

1 24. The system as claimed in claim 22 wherein the correction data 

2 includes linear correction data, 

1 25. The system as claimed in claim 2i wherein the 

2 thermodynamic parameters include DNA thermodynamic parameters. 

1 26. The system as claimed in claim 25 wherein the DNA 

2 thermodynamic parameters include dangling end parameters. 

1 27. The system as claimed in claim 25 wherein the DNA 

2 thermodynamic parameters include coaxial stacking parameters. 

1 28. The system as claimed in claim 25 wherein the DNA 

2 thermodynamic parameters include terminal mismatch parameters. 

1 29. The system as claimed in claim 21 wherein the 

2 thermodynamic parameters include RNA thermodynamic parameters. 

1 30. The system as claimed in claim 21 wherein the 

2 thermodynamic parameters include hybrid DNA/RN A thermodynamic parameters. 
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1 3L The system as claimed in claim 21 wherein the 

2 thermodynamic parameters include DNA loop thermodynamic parameters. 

1 32. The system as claimed in claim 21 wherein the hybridization 

2 information represents top and bottom strand sequences which form a duplex and 

3 wherein the hybridization thermodynamics are calculated for the duplex. 

1 33 . The system as claimed in claim 21 wherein the hybridization 

2 information represents at least a section of a target and a length of at least one 

3 primer or probe complimentary to the target. 

1 34. The system as claimed in claim 33 wherein the hybridization 

2 . thermodynamics are calculated for a plurality of primers or probes complimentary 

3 to the target. 

1 35. The system as claimed in claim 21 wherein the hybridization 

2 information represents at least a section of a target and a primer or probe. 

1 36. The system as claimed in claim 35 wherein a length of the 

2 target is longer than a length of the primer or probe and wherein the hybridization 

3 thermodynamics are calculated for a best target/primer or target/probe complex and 

4 for competitive mismatch complexes. 

1 37. The system as claimed in claim 34 wherein hybridization 

2 information represents at least a section of a target and a primer or probe and 

3 wherein a length of a target is longer than the length of the primer or probe and 

4 wherein the hybridization thermodynamics are calculated for a best target/primer 

5 or target/probe complex and for competitive target/primer or target/probe 

6 complexes. 

1 38. The system as claimed in claim 22 further comprising means 

2 for calculating concentration of each species in a solution at a plurality of 

3 temperatures. 
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1 39. The system as claimed in claim 38 wherein hybridization 

2 information also represents a primer or probe and wherein the length of the target 

3 is longer than a length of the primer or probe and wherein the hybridization 

4 thermodynamics are calculated for a best target/primer or target/probe complex and 

5 for competitive mismatch complexes and wherein the system further comprises 

6 means for calculating concentration of every species in a solution at a plurality of 

7 temperatures, 

1 40. The system as claimed in claim 39 wherein the hybridization 

2 thermodynamics are calculated for at least two best target/primer or target/probe 

3 complexes and for their corresponding competitive mismatch complexes and 

4 wherein the system further comprises means for correcting for any interactions 

5 between the at least two best target/primer or target/probe complexes and their 

6 components. 

1 41. A computer-readable storage medium having stored therein 

2 a database of thermodynamics parameters and a computer program which executes 

3 the steps of: 

4 receiving hybridization information which represents at least one 

5 sequence; 

6 receiving correction data; 

7 receiving a first set of data which represents hybridization conditions; 

8 and 

9 calculating hybridization thermodynamics including net hybridization 

10 thermodynamics based on the hybridization information, the thermodynamic 

11 parameters, the correction data and the first set of data. 

1 42. The storage medium as claimed in claim 41 wherein the 

2 hybridization thermodynamics of individual single stranded, bimolecular and higher 

3 order complexes are statistically weighted in a numerical process and the 

4 equilibrium concentration of each species is output. 
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1 43. The storage medium as claimed in claim 42 wherein the 

2 correction data includes folding correction data. 

1 44. The storage medium as claimed in claim 42 wherein the 

2 correction data includes linear correction data. 

1 45. The storage medium as claimed in claim 41 wherein the 

2 thermodynamic parameters include DNA thermodynamic parameters. 

1 46. The storage medium as claimed in claim 45 wherein the DNA 

2 thermodynamic parameters include dangling end parameters. 

1 . 47. The storage medium as claimed in claim 45 wherein the DNA 

2 thermodynamic parameters include coaxial stacking parameters. 

1 48. The storage medium as claimed in claim 41 wherein the DNA 

2 thermodynamic parameters include terminal mismatch parameters. 

1 49. The storage medium as claimed in claim 41 wherein the 

2 thermodynamic parameters include RNA thermodynamic parameters. 

1 50. The storage medium as claimed in claim 41 wherein the 

2 thermodynamic parameters include hybrid DNA/RNA thermodynamic parameters. 

1 51. The storage medium as claimed in claim 41 wherein the 

2 thermodynamic parameters include DNA loop thermodynamic parameters. 

1 52. The storage medium as claimed in claim 41 wherein the 

2 hybridization information represents top and bottom strand sequences which fonn 

3 a duplex and wherein the hybridization thermodynamics are calculated for the 

4 duplex. 
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1 53. The storage medium as claimed in claim 41 wherein the 

2 hybridization information represents at least a section of a target and a length of at 

3 least one primer or probe complimentary to the target. 

1 54. The storage medium as claimed in claim 53 wherein the 

2 hybridization thermodynamics are calculated for a plurality of primers or probes 

3 complimentary to the target. 

1 55. The storage medium as claimed in claim 41 wherein the 

2 hybridization information represents at least a section of a target and a primer or 

3 probe. 

1 56. The storage medium as claimed in claim 55 wherein a length 

2 of the target is longer than a length of the primer or probe and wherein the 

3 hybridization thermodynamics are calculated for a best target/primer or target/probe 

4 complex and for competitive mismatch complexes. 

1 57. The storage medium as claimed in claim 54 wherein 

2 hybridization information represents at least a section of a target and a primer or 

3 probe and wherein a length of a target is longer than the length of the primer or 

4 probe and wherein the hybridization thermodynamics are calculated for a best 

5 target/primer or target/probe complex and for" competitive target/primer or 

6 target/probe complexes. 

1 58. The storage medium as claimed in claim 42 wherein the 

2 program further executes the step of calculating concentration of each species in a 

3 solution at a plurality of temperatures. 

1 59. The storage medium as claimed in claim 58 wherein 

2 hybridization information also represents a primer or probe and wherein the length 

3 of the target is longer than a length of the primer or probe and wherein the 

4 hybridization thermodynamics are calculated for a best target/primer or target/probe 

5 complex and for competitive mismatch complexes and wherein the program 
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6 executes the step of calculating concentration of every species in a solution at a 

7 plurality of temperatures. 

1 60. The storage medium as claimed in claim 59 wherein the 

2 hybridization thermodynamics are calculated for at least two best target/primer or 

3 target/probe complexes and for their corresponding competitive mismatch 

4 complexes and wherein the program executes the step of correcting for any 

5 interactions between the at least two best target/primer or target/probe complexes 

6 and their components. 
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Figure 1 
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Hybridization Information 

Top Strand/Bottom Strand Types 
DNA/DNA {▼] 

Top Strand Sequence 5-3' 

cccaaaaaaaaaaaccg 



Module 1 



Bottom Strand Sequence 5'-3' p] 
J7j | Use Complement 



*ggtttttttttttgg* 



J3 



JB 



Hybridization Conditions 

: User defined values for (Na+] and [Mg2+j J7j 
[Monovalent cation] ; o.ios | 



[Mg 2+ ] 

Hybridization 
Temperature 

[Top Strand ] 
[Bottom Strand] 



J 



5e-8 



3e-7 



mol/L 
mol/L 

°C 

J mol/L 
J mol/L 



Corrections 

Linear Correction for Micro Chips 
(AG 0 37 (microchip) = ax AG° 37 (solution) + b) 

a= XZJ b = J ) 

Top Strand Folding Correction 

AG °37= 3L1 J kcal/mol AH°= -37.8 \ kcal/m 
Bottom Strand Folding Correction 

AG °37= 0 1 kcaVmol AH°= J> j kcal/m 



| Predict Thermodynamics j Clear Input 



Figure 2a 
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Duplex sequence 

5 / -CCCAAAAAAAAAAACCG-3 / 
3 / -*GGTTTTTTTTTTTGG*-5 / 

Experimental conditions 



Hybridization type = DNA/DNA 
[Top strand] = 0.5E-07 mol/L 
[Bottom strand] = 0.3E-06 mol/L 

Hybridization temperature = 37.0 °C 

Corrections 



Top strand folding: 
AH° = -37.8 kcal/mol 
AG° 37 0 = -2.10kcaymol 



f Module 1 . 



Thermodynamic predictions 



In 1.000 MNaCl: 

A H° = -119.3 kcal/mol 

AS 0 = -335.8 eu 

AG° 3? 0 = -15.14 kcal/mol 

T M = 52.9°C 

In 0.1050 M NaCl and 0.0000 M 
MgC12: 

A H° = -119.3 kcal/mol 
A S° = -348.3 eu 
AG° 370 = -11.29 kcal/mol 

T M = 42.2°C 

The net hybridization thermodynamics 
is: 

AG° 370 = -8.74 kcal/mol 
T M = 34.9 °C 



Note: 

The net hybridization temperature is the temperature at which the concentration of duplex equals half the 
maximum possible concentration of duplex. 

The net free energy is calculated from the net equilibrium constant at the given temperature: 
Knet=pDuplex]/((Ct-[Duplex])*(Cb-[Duplex])), where [Duplex] is the concentration of duplex, Ct is the 
initial concentration of top strand, Cb is the initial concentration of bottom strand. 



Figure 2b 
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Hybridization Information 

Target/Primer Types 
DNA/DNA [▼) 

Target Sequence So* 



Module 2 



Primer length 15 1 
Number of best primers to be displayed J* j 



Fj 



Hybridization Conditions 

'.A ...«.AAAi AAA V .A.. . ■ . >. ._• . .aw 

' User defined values for [Na+] and [Mg2+] |t| 



[Monovalent cation] 


: 1 


i 


mol/L 


[Mg 2 *] 




i 


moI/L 


Hybridization 








Temperature 




°C 


[Target] 






^| mol/L 


[Primer] 






] mol/L 



Corrections 

Linear Correction for Micro Qiips 
(AG° 37 (microchip) = a X AG° 37 (solution) + b) 



a = 



J 



b = 



Figure 3a 
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Module 2 



Target sequence 



5UCCGTTTGTA GTCCGTACGA CAC AT AACGG TGCATTC 



[Na + ] = 1.0000 mol/L 
\M^ + ] = 0.0000 mol/L 

The 2 best primers of length 15 are: 



From position 28 to 42: 5 -GGTTGCAATGCACCG -3' 

AH° = -132.0kcal/mol AS° = -355.6 eu AG 0 37 0 = -21.71 kcaVmol T M = 
70.2 °C 

From position 35 to 49: S'-GCAGCATGGTTGCAA -3 ; 

AH° = -124.8kcal/mol AS 0 = -336.5 eu AG 0 37 0 = -20.42 kcal/mol T M = 
68.4 °C 



Experimental conditions 



Corrections 



Hybridization type = DNA/DNA 
[Top strand] = 0.1E-05 mol/L 
[Bottom strand] = 0.1E-05 mol/L 

Hybridization temperature = 37.0 °C 



No corrections 



Figure 3b 
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Module 3 



Hybridization Information 

Target/Primer type 
PNA/PNA [^| . 

Target 5'-3' Primer Sequence 5'-3' , [7] 

• acgcttgaatgcagttaatgcc £] tgaatgcagt ' ' □ 



Minimum percent stability of alternative binding, sites 
compared to the most stable binding site 

Number of base pairs required to compute the solution jj[ZLl ' I 



L_ B 



Hybridization Conditions 

■ User defined lvalues for {Na+j and (Mg2+j Jt| 



[Monovalent cation] Hl 
[Mg 2+ ] £o 

Hybridization 
Temperature 

[Target] 

[Primer] 



2 



J 



37Vcf 



le-6 



mol/L 
mol/L 

°C 
J mol/L 



Corrections 

linear Correction for Micro Chips 
(AG 0 37 (microchip) = a x AG° 37 (solution>+ b> 

a = j. I t>= j i 



j Submit"" | Clear Input 



Figure 4a 
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Module 3 

Target sequence 



5'-ACGCTTGAAT GCAGTTA ATG CC-3' 

Primer sequence 



I'-TGACGTAAGTS 1 

Experimental conditions Corrections 



Hybridization type = DNA/DNA No corrections 

[Top strand] = 0.1E-05 mol/L 
[Bottom strand] = 0.1E-05 mol/L 

Hybridization temperature = 37.0 °C 

[Na + ] = 1.0000 mol/L 

[Mg 2+ ]= 0.0000 mol/L 

Number of base pairs required to compute the 

solution = 5 

Best primer site 



from target position 8 to position 17 

GAATGCAGTTAA 
TGACGTAAGT 



A H° = -26.2 kcal/mol AS° = -70.7 eu AG° 37 0 = -4.28 kcal/mol 
T M = -9.9 °C 

Figure 4b 
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Module 5 



Hybridization Information 

Target/Primer Types 
DNA/DNA [rj 

Target Sequence 5'-3' 



B 



Find best primer in sequence section ranging from 
nucleotide number: 



^1 



J 



Primer length 
p| Number of best primer 



to 

. L IU l 



10 



J 



Percent stabilitty of alternative binding sites compared to the most stable binding site so ) 
Number of base pairs required to compute the solution 7 J 



Hybridization Conditions 

• User defined values for [Na+j and [Mg2+] frj 
[Monovalent cation] fi ~ jj mol/L 

tMg 2 ^ ro VV 1 mol/L 
Hybridization :•-»«•■■ ^ o 
Temperature ' '3 C 

[Target] rje-6 ' ' j mol/L 

[Primer] ' i«-6 j mo l/L 

j Submit 



Corrections 

Linear Correction for Micro Chips 
(AG° 37 (microchip) = a X AG° 37 (solution) + 

±—1 0 I 



a = 



j Clear Input 



Figure 5a 
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Target sequence 

• ModuleS 



5'- 

AGGTCCATGCTTTGGAACAGCTACTTGAACCGATCATGGACACTGACGGATAAC( 

-y 

Experimental conditions Corrections 



Hybridization type = DNA/DNA 

[Top strand] = 0.1E-05 moI/L No corrections 

[Bottom strand] = 0.1E-05 mol/L 

Hybridization temperature = 37.0 °C 

[Na+] = 1.0000 mol/L 

[Mg 2+ ]= 0.0000 mol/L 

Number of base pairs required to compute the 

solution = 7 

Best primer search area from position 1 to 
position 60 



Best primer # 1: 



from target position 35 to 49 
5 1 -TCATGGACACTGACGGA- 3 1 
3 1 -GTACCTGTGACTGCC - 5 / 

A H° = -123.1 kcal/mol AS 0 = -331.5 eu 
68.4 °C 



AG 0 yj 0 = -20.27 kcaymol T M = 



Best primer #2: 



from target position 18 to 32 
5 / -ACAGCTACTTGAACCGA- 3 f 
3 f -GTCGATGAACTTGGC - 5 1 

•A H° = -125.0 kcal/mol AS 0 = -339.6 eu AG 0 
66.1 °C 



370 = -19.67 kcal/mol T M 
Figure 5b 
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initial guess: 
C s$p =C p To,a, /2,C sst = C t T<> «» , /2 
temperature T=0 °C 



Iterativery solve for C sst and C ssp 
(Primer and Target conservation equations) 



NO 




YES 




YES 



T=T+1 



NO 



Figure 6 



Multiplex PCR Design 



PI 



> Gene I 



P3 ^ tt Target DNA 

r Gene II 



JCI1C 1 C. i .... I L 

I I 1 M i i 1 M l I I I 1 II I M 1 3 



P2 



P4 



Figure 7 
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Prediction of Molecular Beacon 
Hybridization 



Duplex Formation 

CCCAAAAAXAAAAACCG 
+ 

GGTTTTTYTTTTTGG 



X-Y 



CCCAAAAAXAAAAACCG 
GGTTTTTYTTTTTGG 



Beacon Folding 



Random Coil Beacon 



AG° 37 =-2.1 kcal/mol 



T M = 55.2°C 



Net Hybridization 

Hairpin Beacon 
+ 

Target 



A-T -13.24 47.9 

A-A -9.94 39.1 

A-C -9.03 36.0 

A-G -10.36 40.4 

A X a 
A A A A 

A A 

A A 

A A 

C C C C 
C-G 

T-A 

C-G 

G-C 

C-G 

•) 

Hairpin Beacon 



©-CGTCCCCAAAAAAAAAAACCGACG-C 

~ / | \ 3 ' -GGTTTTTXTTTTTGG-5 ' 



Target 



X-Y 


AGV 


(Effective) 


T M (Effective) 




Exp. 


Pred. 


Exp. Pred. 


A-T 


-10.49 


-10.69 


42 42.4 


A-A 


-6.66 


-7.39 


27 26.8 


A-C 


-6.72 


-6.48 


23 21.1 


A-G 


-7.62 


-7.81 


28 29.5 



3 GG I 1 1 1 1 1 I I I I I GG 5 ' 
3'GGTTTTTAI I I I IGG 5 ' 
3 GG I I I I I C I I I I I GG 5 ' 
3 'GGI I I I IGI I I I IGG 5 ' 



0.105 MNaCl 0.001 MMgCl 2 [beacon] = 5xl0 8 M [target] = 3x10"^ 
Bonnet et al. (1999), Proc. Nat. Acad. Sci. USA 96, 6171-6176 



Figure 8 
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Simulation of Molecular Beacon Hybridization 

1.2 1 : 1 
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Match vs. Mismatch Hybridization 

Match Site Mismatch Site 
1 I I I ■ p—i 1— I Target DNA 

MIT 11x11 

Probe Probe MM 

Equilibria 

Match Hybrid T + P + 
Mismatch Hybrid T + P< 
Double Hybrid T + 2P ^ 

• Given Cr^ftotal], C Probe [total]and the 3 equilibrium constants above, 
it is trivial to solve for the concentrations of all species 

• Since AG° 37 and AH 0 are known, calculate K's at all temperatures 

• Simulate hybridization at all temperatures - optimize specificity 

• More complex model would also include single-strand folding equilibria 

FigurelO 



± T-P 

"p.pMM 

dt Y_p_pMM 



K M = 



K DH~ 



[T-P] 

[T] [P] 

[T-P MM ] 
[T][P] 

[ T-P-P MM ] 
[T][P] 2 
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Match vs. Mismatch Hybridization Simulation 



c 
o 

2 
*■• 
c 

V 

u 
c 
o 
o 

0) 
N 

"5 
E 



1.00 
0.90 
0.80 
0.70 
0.60 
0.50 
0.40 
0.30 
0.20 
0.10 
0.00 



9 




♦ Fraction Target in R.C. 
■ Fraction Probe in R.C. 
£ Fraction Match Hybrid 
O Fraction Mismatch Hybrid 
▲ Fraction Double bound 



80 



100 



Temperature ( C) 



Figurell 
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